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ABSTRACT 


The present thesis reports studies on two aspect 
discharge measurement devices, viz. (i) the selection 
flow measurement devices and design and (ii) the measure 

• rpv> p g Y "* 0 

ment of discharge in submerged flow condition. 

reported in part A and part B of the thesis respecti Y 

The first study reported in part A concerns with rele 
vant devices by considering the constraints and attribute 

« _ i s 

of twenty six open channel flow measuring devices, 
is achieved through a database program SELECT, written 
FORTRAN IV. For the design of flow measuring devices a 
program DESIGN is developed. This program uses trial and 

error procedure to develop the weir geometry. 

. t-v no on an 

The second study# reported in part B# i 

experimental investigation of the submerged flow m six 
kinds of weir shapes. In the generalised form of Villem 
equation the variation of Kg and m are obtained for 
shapes tested. For convenience of discharge prediction two 

zones of submergence are suggested as follows*. 

,, j 0 c. f a 50 percent 

i) Zone I. Where the submergence is o 

. „ A n fo 95 percent. 

ii) Zone II. Where the submergence is 

The values of K s and m are obtained for all th 
weirs# in these t Wo zones. Equations for the 
prediction of submerged flow discharge with an 



accuracy of 97 percent in Zone I and 95 percent in 
Zone II are obtained. The modular limit at five 
percent level (M^) and ten percent level (M^q) 
are also obtained for each of the weirs. These 
relationships enable submerged flow mode to be 
used for discharge measurement in sharp crested 
weirs at an accuracy of about 95%. 
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PART- A 

The selection of flow measurement device 
and its design 



CHAPTER A.l 
INTRODUCTION 


A. 1.1 Introduction 

The efficient measurement of discharge in channels is 
very important for the proper utilization and management of 
available water resources for irrigation# water supply and 
waste water disposal etc./ so that we can fulfil the obje- 
ctive of available water resources according the demand of 
users. Also accurate flow measurement is a must for proper 
charge to be fixed for quantity (volume) of water used or 
effluent discharged into a stream. 

The measurement of discharge in an open channel can be 
done in two different modes of flow viz. free flow and 
submerged flow. Both the flow regimes have their own attri- 
butes. If the discharge measurement is done in free flow 
condition# only the upstream head has to be measured. On 
the contrary# for submerged flow discharge measurement, 
the measurement of both the upstream and downstream heads above 
the crest level are necessary. 

In the present study two aspects of discharge measuring 
devices are considered. These are 

(i) the selection of flow measurement device and its 
design 

and (ii) the flow measurement under submerged flow 


condition 
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These are described in two parts, viz. in part A and Part B 
respectively. 

In this chapter, the description of the various flow 
measuring devices commonly used in practice are given. 

A. 1.2 Principle of flow measuring device 

Flow Measuring Devices: The commonly used open channel flow 
measuring devices are considered under four categories as, 

1. broad crested weirs 

2. flumes 

3. sharp crested weirs 
and 4. orifices. 

The basic principle of discharge measurement for broad cres- 
ted weirs and flumes are based on the critical depth theory 
and their relationship with total energy. To achieve the 
critical depth in an open channel flow at a geometrically 
specified constriction is built with the help of broad cres- 
ted weir and flumes at the crest of broad crested weir and 
in the throat of flume respectively where sufficient fall 
is available. Since there is fixed relationship between 
the specific energy and the critical depth, a single depth 
in free flow enables the discharge to be estimated in these 
devices. 

In the case of sharp crested weirs and the orifices, 
the water approaching them with relative low velocity acce- 
lerates to high velocity as it issued out from the crest of 
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a weir or through an orifice. Based on total specific energy# 
a relationship for the discharge in term of the head over a 

weir or an orifice is established. 

A.1.3 Types of flow measuring device 

The commonly used measuring devices in both the cate- 
gories can be classified into twenty six different catego- 
ries as below based on the shape and operating characteristic. 

1. Broad crested weirs and sills 

1.1 Round nose horizontal broad crested weir 

1.2 Romjin movable weir (Ref . l) 

1.3 Triangular broad crested weir 

1.4 Broad crested rectangular profile weir 

1.5 Faiyum weir (Ref. l) 

1.6 Weir sill with rectangular control section 

1.7 V-notch weir sill 

1.8 Triangular profile two-dimensional weir 

1.9 Rectangular profile flat-Vee weir 

1.10 Butcher's movable standing wave weir 

1.11 WES-standard (ogee) spill way (Ref. l) 

1.12 Cylindrical crested weir (Ref. 1). 

2. Flumes 

2.1 Long throated flume 

2.2 Parshal flume 

2.3 H— flume. 

2.4 Throatless flume 
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3. Sharp Crested Weir 

3.1 Rectangular sharp crested weir 

3.2 V-notch sharp crested weir 

3.3 Circular sharp crested weir 

3.4 Parabolic sharp crested weir 

3.5 Proportional weir. 

4. Orifices and Sluice Gate 

4.1 Rectangular sharp-edged orifice 

4.2 Radial or tainter gate 

4.3 Crump-de Gruiter adjustable orifices (Ref. 1) 

4.4 Meter gate 

4.5 Neyrpic gate (Ref. 1). 

Each has a specific range of various characteristics in which 
it is operational or efficient. 



vunri iiK Z 
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CHARACTERISTICS AND PROBLEMS 


A. 2.1 Characteristics : The following are the characteristics 

of a flow measuring device* 

(1) The shape of control section perpendicular to flow 
direction. 

(2) Possible functions of a structure. 

(3) The minimum and maximum values of upstream head. 

(4) The minimum and maximum values of downstream head. 

(5) The minimum crest height above the approach channel 
bed. 

(6) The minimum size of control width. 

(7) The range of discharge measurement. 

(8) The sediment discharge capacity. 

(9) Floating and suspended debris passing capacity. 

(10) Modular limit. 

(xi) Flexibility. 

Taking them in detail. 

(1 ) The shape of control section perpendicular to flow 
direction* 

It depends on head and discharge relationship and 

the maximum and minimum discharge. If variation in the 

Q 

maximum and minimum discharge (i.e. Y = "— ■ ■■ ) are high 

Hnin 

the triangular shape of control section is suitable. If 
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3 

maximum discharge is less than 0.4 m /sec and the sediment 
and trash are not necessary to be passed then the sharp 
crested weirs are suitable. Further the shape of sharp 
crested weir depends on the head discharge relationship. 

3 

For high discharge above 6 m /sec/m the shape is controlled 
by various practical and constructional reasons. So Ogee 
spillway or cylindrical crested weirs are suitable for 

3 

Q = 6 m /s/m. 

(2) The function of measuring devices * 

The measuring devices may have one or more function (s). 
They are 

(a) measurement of discharge without discharge 
regulation 

(b) measurement of discharge with discharge 
regulation 

and (c) the volumetric measurement of discharge. 

If the area of control section can not be changed, the 
structure can only be used to measure discharge without 
regulation. On the contrary# if the area of control 
section can be changed with the help of a movable gate 
arrangement or movable crest level# the measurement of 
discharge can be done with discharge regulation. The 
propeller meter (Ref. 1) and the dethridge meter (Ref.l) 
can measure a flow rate and totalise the volume in cubic 
meters. In this condition discharge can be regulated 
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with the help of separate movable gate. So the functional 
requirement of flow measurement with and without discharge 
regulation is an important parameter in the selection of 
relevant discharge measuring devices. 

(3 ) The minimum and the maximum upstream head* 

The minimum upstream head, m j_ n ' and the maximum 
upstream head H 1 are very important parameters for the 

validity of head discharge equation used for any particular 
measuring devices. The value of discharge coefficient, C d 
and velocity coefficient, C v are also effected by the upstream 
head. Since each type of measuring device has his own fixed 
range of head within, which it is operational and efficient. 
The upstream maximum head is also limited because of various 
practical and constructional reasonsc 

(4) The maximum head of downstream: 

The maximum downstream head may produce submerged flow 
condition through a discharge measuring device. The sub- 
merged flow behaviour of different measuring devices are 
different. Some-times the downstream head of measuring 
device may increase the water level of the upstream and water 
may over flow through the embankment of channel. It also 
can produce error in discharge estimation which is mentioned 
in modular limit. 
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(5) The minimum crest height above the approach channel ! 

The constriction is applied with the help of rising 
the crest level above the bed of channel in a measuring 
device to obtain required flow situation for their head 
discharge relationship. Since the required flow situation 
for flow measurement varies with measuring devices, so 
the minimum value of crest height above the bottom channel 
are different for different measuring devices. 

For the water supply channel (in irrigation) and for 
feeding water from main canal to branch canal the minimum 
water level is also maintained by rising the crest level. 

To fulfil these requirement, selection of proper value 
of P and selection of corresponding efficient measuring 
device is necessary. 

(6) Minimum size of control width; 

The width of control section controls the upstream 

3 

head and the discharge in meter /sec/meter. It also 
affect the critical depth in the constriction zone of meas- 
uring devices and the water level in the upstream zone. 

So for a particular situation of open channel flow the 
judgement of suitable control width with relevant measur- 
ing devices are necessary for their efficient operational 
behaviour. 



(7) The range of discharge measurement: 

The minimum and maximum discharges are very important 
characteristics for a measuring device for their efficient 
operational characteristics. For the maximum discharge 
the selection also depends on various practical and cons- 
tructional reasons. Since the minimum discharge effects 
the minimum head. For each type of measuring device, 
there is a limiting value of head up to which their head- 
discharge equation is applicable. For examples 


In case of broad crested rectangular profile weir 


H 1 max 

0.6 



H 0 

2 max 

0.06 

For P = 

1.0m 

Sim ax ” 

1.0 

L = 

0.5 m 

C dmin ~ 

0.84 




where L = length of broad crested weir along the flow 
direction. 

Then, 

max _ Y 
Q 1 min 

Y 


"dm ax 
-* 

"dmin 


H, 1. 5 

x C-J-ESX, 


H 


1 min 


1.0 

0.84 


x 


(■ , Q . a.§~ ) 

V 0. 06' 


1.5 


= 37 


(8) The sediment discharge capacity : 

The weir, sill, flumes and tainter gate are efficient 
discharge measuring structures where the sediments discharge 


are necessary 
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(9) Floating and suspended debris passing capacity: 

The sharp crested weir and orifices are easily clogged* 
So they are not recommended if floating debris has to be 
passed. The weir with rounded nose and flumes are efficient 
for passing the floating and suspended debris. 

(10) Modular limit : 

The modular limit of a given measuring device is 
defined as the value of submergence ratio Ir^/H^ at which 
discharge deviates by one percent from the discharge calcu- 
lated by free flow head discharge relationship. So for a 
given situation of flow through a channel the ratio of 
maximum downstream head and upstream head should be consi- 
dered as important quantity in the selection of discharge 
measuring devices with reasonable accuracy. A brief review 
of modular limit is presented in section (B.3.4). 

(1 1 ) FI exibility t 

The flexibility is an important characteristic for 
flow measuring structures at the junction of main channel 
and off take channel where flow regulation is necessary to 
distribute the available discharge between off take channel 
and main channel according the requirement of available 
water in a channel. The flexibility of an off take channel 
measuring device is defined as the ratio of the rate of 
change of discharge in the off take channel to the rate of 
change in discharge of the main channel. The flexibility 
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where Q m = discharge passing through main 

channel « 

The value of the flexibility may be less than one, equal to 
one and greater than one. Depending on the value of flexi- 
bility a measuring device at the junction of main channel 
and off take channel, the measuring device may be sub- 
proportional, proportional and hyperproportional respectively » 
By selecting a proper value of the height of crest level 
and proper category of discharge measuring device the required 
flexibility of off take measuring device can be obtained 
and vice-versa. The following example illustrates this 


point : 
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Example 1 * An irrigation main channel contains an H— flume 
for discharge measurement. Select a suitable branch channel 
measuring device having fixed crest at the level of the 
floor of the H-flume. In this problem the given data ares 

= discharge ratio for main channel = = 2 


“•m 




xn 


°m 


Y = discharge ratio for off take channel = 

° Sm 


ax 


Solution : 


As 


xn 


8 

7 


dQ m /Q , 


max 


- Q. 


min 


m 


dQ o /D o 


Q 


max 


max 


°min 


1 

2 

1 

8 


flexibility = F = 


d °o /Q o 


• 3^75 


m 


= I X I 

8 I 


1 

I 


n x H 

„ o m 

K SB . . i i — # 

n x H 

m o 


since the crest heights are same (i^ = H Q ) 


As for H-flume n^ = 2, 


F 



1 

4 ' 


n 

o 


1 

2 * 


Hence# a rectangular sharp edge orifice is recommended. 

Table A. 2.1 gives a summary of the characteristics of 
the 26 devices used in the present study. 



Summary of the Characteristics 
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A. 2.2 Factors Affecting Selection of a Device : 

The different characteristics of a flow measuring 
device (as in section A. 2.1) are interrelated with one 
another and these characteristics have their own limita- 
tion. So the selection of an ideal flow measuring device 
for a particular flow situation to fulfill all the funct- 
ional requirement is difficult. The relevant few structures 
can be selected by comparing the function requirement with 
the characteristics of flow measuring devices. The major 
functional requirements are as follows: 

(i) discharge measurement only 

(ii) discharge measurement with regulation 

(iii) the required head loss 

(iv) required upstream and downstream water level 

(v) the accuracy in discharge measurement 

(vi) required flexibility to fulfill the demand of 
water according their necessity 

(vii) to pass sediment or not 
and (viii) to pass trash and floating material. 

These requirements can be fulfilled on the basis of the 
eleven important characteristics of flow measuring devices 
listed in Section A.2.1. As many as 26 different commonly 
used measuring devices are considered inthis study. Each 
one of them has its own attributes and range of characteristics 



A very large number of factors influence the selection of 
proper measuring device. In a practical situation an 
engineer facing problem of selection of a proper meas- 
uring structure/ should have proper knowledge of attri- 
butes of each type of structure so that he can compare 
with existing situation and functional requirement. But 
due to lack of requisite knowledge and systematic analysi 
of attributes# the probability of a choice of improper or 
inadequate flow measuring device is indeed large. So 
the present study deals a method to use the computer to 
aid in the selection of a set a feasible devices based on 
a database program for a particular problem. The details 
of the program along with typical examples are presented 
in Section A. 3. 



CHAPTER A. 3 


SELECTION AND DESIGN 

A. 3*1 Program: SELECT 

The program SELECT is based on database which follow 
a well defined algorithm that specifies explicitly how to 
find out variables for any given input variables. The 
output variables are sorted out step by step on the basis of 
a given set of input variables. The flow chart showing the 
various components and functions of the program SELECT is 
shown in Fig. A. 3 . 1 . 

The rules of the program is fixed initially on the 
basis of the following requirements: 

(1) Prime function of measuring structure (with or 
without discharge regulation) 

(2) Sediment to be passed or not 

(3) The trash and floating material have to be 
passed or not. 

The following input data are required for the selection 
criteria for discharge measuring devices. 

(1) The maximum discharge = C^ ax in Cumec 

(2) The minimum discharge = Q . in Cumec 

mm 

(3) The maxim-urn depth of water flow = in meters 

(4) The minimum depth of water flow = Y m ^ n in meters 

(5) The width of control section = B in meters 

(6) The crest height above the channel bed = P in meters 

(7) The required maximum head loss i n meters 




FIG. A. 3.1 : Flow chart of program SELECT 

Fig. -A. 3. 1.1 









Section* 
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and (8) The flexibility of a measuring structure at the 
junction of two or more channels* 

The first four data can be obtained directly by field 
investigation of an open channel flow. The data (5), (6) 

and (7) depend on the boundary required to fulfill other 
objectives. The data (8) is required only when the flow 
takes place at the junction of two open channels. On the 
basis of the above data, following parameters are calculated s 


(i) 

Discharge ratio 

= Y 

max 

^in 

(ii) 

Maximum upstream 

head 

= H. ^ ^ in meters 
1 max 

(iii) 

Minimum upstream 

head 

= H. 

1 mm 


(iv) Maximum downstream head = H 0 

2 max 

and (v) The exponent of head discharge equation = n. 

On the basis of these input values, output values are sorted 
out step by step and the desired relevant structure (s) is/are 
obtained. 

The program SELECT is given in Appendix AA-1. This 
program is written in FORTRAN IV. It is simple, efficient 
and interactive in nature. It can be handled on a personal 
computer (PC) with small memory. 

Typical solution obtained by using program along the 
respective input is shown in Append ix-AA-2« 

It is seen that the program gives a set of possible 
devices that could be used in a circumstances. The decision 
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to adopt any one of this set depend upon the case consider- 
ation in fabrication and maintenance, environmental and 
factors peculiar to the fabrication. 


A. 3.2 Program DESIGN 

After selection of the relevant measuring devices 
their design can be done on the basis of input data by 
trial and error procedure. Each type of measuring device 
will have definite head discharge relationship. The actual 
value of the relevant coefficient of discharge, or the 
coefficient of velocity C v# depends on the geometry, the 
shape and size of a given device and its installation 
(Refs. 1, 3) . 

Initially trial value of C d and C v are assumed and 
the corresponding value of the heads over the crest are 
found out. Again the next trial values of and C v are 
found out and their corresponding heads are calculated. 

The process is continued upto desired accuracy of heads 
above crest level. On the basis of head over the crest 
the actual required value of crest height is calculated. 

All the other standard values were calculated, for example 
(Ref. 3). In the case of Ogee spillway having vertical 
face of upstream spillway the values of the co-ordinates 
were calculated by equation. 


5.31 


*1 Q C 

0.50 (■■ — "T — ) * for downstream profile 


(A. 3.1) 


5.31 ~ 0.724 ( 5>31 + 0.27) 


0.432 


/ x 


0.625 
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The computer program DESIGN is developed to provide 
the basic design of the following type of weirs: 

(1) Round nose broad crested weir 

(2) Romjin movable measuring/regulating weir 

(3) Triangular broad crested weir 

(4) Broad crested rectangular profile weir 

(5) WES-standard spillway 

(6) Linear proportional weir (Ref. 3) 

(7) Quadratic weir (Ref. 3) 

(8) Logarithmic weir (Ref. 3) 

(9) Exponential weir (Ref. 2) 

(10) New baseless weir (Ref. 2) 


The relevant constraints and the coefficient of discharge 
of each type of the above weir shape is given in Table A. 3.1. 

The program is written in FORTRAN IV and given in 
Appendix AA-3. A few selected design developed by this 
program along with their respective inputs are shown in 
Appendix AA-4 . 

Other Devices: The design of different types of propor- 
tional weirs viz. quadratic weir, exponential weir, logari- 
thmic weir and new baseless weir were considered for the 
design which are not considered in the selection of discharge 
measuring devices. The consideration of these weirs were 
done to the complex nature of their profile which is 
difficult to and time consuming design. 
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For example the profile equation for exponential 
weir is 


= b fi - 


| tan ' 1 ® 1 / 2 

it a 


? + 


2K 3/ ^ 2 e Y 


( 2 g) 


1/2 


* e~ z dz 


o 

(A c 3 .1) 


where 


a - K e- Va 

(A. 3. 2) 

K = j (2b) c d (2g ) X / 2 a 3 / 2 

(A .3 .3) 

^ =!a 

(A.3 .4) 


where \ is datum parameter. 


(ii) The logarithmic weir profile 
x = b {1 - r tan 1 ( Y ) 1//2 ? - — 


(y/a) 


1/2 


71 


where. 


( 6 n + y/a) 


<^n 

(<Tn 

1 r, 

Irfn + y/a) 3 ^ 3 

-LJ.I 

(<fn 

Q = K In ( 6 + 

-5-) 


na 

2 b(2g) 1/2 C d | 

a 3 / 2 : 

2b(2g) 1/2 C d a 172 

= 


(tfn + y/a) + (y/a) 


K 

/na 


1/2 

T 172 

(A. 3. 5) 
(A„3 . 6 ) 

(A. 3. 7) 
(A. 3. 8 ) 


where b is datum parameter, n is an arbitrary number, 

a and b are height and half of the width of base weir 
respectively.. 
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Equations (A.3:.7) and (A. 3. 8) may be used to solve 
the values of any two of the parameters such as a, C and b 
for an assumed value of n. The point of zero width can 
be shifted by choosing arbitrary higher value of n. So 
the design of logarithmic weir for maximum and minimum 
discharges measurement and base width, based on trial and 
error procedure. These weirs have advantage due to their 
small relative error. 



31 


CHAPTER A— 4 

CONCLUSIONS AND RECOMMENDATIONS 

A. 4.1 Conclusions 

(1) A data based program (SELECT) has been developed 
to identify a set of possible measuring devices for a speci- 
fied set of input data. This program scans 26 devices each 
having upto ten attributes. 

The program SELECT will be of great help to field 
engineers in identifying an appropriate flow measuring 
devices in a given open channel flow situation. 

(2) A program (DESIGN) has been developed to 
provide the relevant design parameters along with relevant 
weir shape profile for ten weir shapes. This program is 
particularly helpful in the design of different types of 
proportional weirs viz. exponential weir* new-baseless weir/ 
logarithmic weir etc. 

A. 4 . 2 Recommendations 

(1) The design of other device shapes listed in 
program SELECT can be done in the similar fashion. The 
"star-base" on HP workstation can be used for schematic 
sketching with detailed design. 
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(2) A combined interactive program can be developed 
for the selection as well as for the design of relevant 
discharge measuring devices. 

(3) An expert system should be also developed on 
the basis of program SELECT. 
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CHAPTER B.l 
INTRODUCTION 


B.1.1 Introduction 

When water is flowing over the sharp crested weir in 
such a fashion that the tailwater level is above the crest 
level of weir then the flow is known as submerged flow. The 
thin plate weirs are extensively used for the measurement 
of small discharges due to their accuracy. Normal ly, these 
are used in free mode only. If the discharge can be meas- 
ured with reasonable degree of accuracy (say within five 
percent of error) in the submerged flow condition also then 
these devices would have larger range of application. Also, 
the following advantages of submerged flow can be obtained. 

(1) The energy loss occurring through free flow weir 
will be as great as the energy loss due to considerable 
length of flow. This loss of energy can be saved consider- 
ably with submerged flow mode of discharge measurement. 

(2) In free flow the floor of the channel will have 
to be increased by larger amount above the channel bed 
there by increasing the silting and seepage loss. But by 
having submerged flow conditions the silting and seepage 
loss can be reduced due to reduction in the crest level 
above the bed of channel. 
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Surprisingly# the submerged flow over sharp crested 
weirs# have not received much attention of research wor- 
kers. Vil lemon te (l) is probably the first to have con- 
ducted a systematic study of submerged flow over weirs. 
Rangaraju (6) reported submerged flow over broad crested 
weirs. However# there is no reported work other than 
that of Villemonte to the best of authors knowledge 
on submerged flow over thin plate weirs. 

A brief review of Villemonte work is given in the next 
section (B.l.2) . 


BJL.2 Villemonte*s Study 

Consider a submerged weir flow as in Fig. B-l# Villemonte 
assumed that the net flow over the weir could be expressed 
as the difference of free flow discharge due to head in 
upstream minus the free flow discharge due to head H 2 in 
the downstream. 

Thus# if Q s = submerged flow# and are free 

flow discharged due to head and H 2 then 



This assumption implies that the head and H 2 do not directly 
affect the flow of water due to head of upstream above 
the crest level of the weir and vice versa. 
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Fig. B.1 Submerged flow weir 



The equation (B..1 . 2) should not be expected to give 
exact quantitative relationship for determining Q g , as 
there are interactions and other perturbing influences 
which are completely ignored. So the equation (B. 1.2) was 
modified by Villemonte as. 


Q. 


= fd - £2) = 

U fl 


K s (1 - 
s Q fl 


Since, 


'f 1 


n i 

C 1 H l ' Q f2 


n, 

c 2 V 


then 


n. 

a, h 0 4 

s 2 2 


Q 

n = K s (1 - 
Q fl S 


n„ 

C 1 H l' 


(Bo 1 . 3) 


(B. 1.4) 


where and C 2 / and n 2 are coefficients and exponents 
that appear in free-flow discharge equation. For a given 
type of weir C ± = C 2 , ^ = n 2 = n, and taking Q fl = Q f 



H 


K s fl - 


m 


Hi 


(B.1.5) 


In this equation K_ and m are to be found by experi- 

s 

ment. Villemonte conducted a large series of experiments 

in 1943 on seven type of weirs with values of n varying from 

1.0 to 3.5. For each weir he assumed K g = 1.0 and mean 

average value of m was calculated by drawing the best fit 

H 2 n Q s 

line on a logarithmic plot of 1 - (rr-) n and — . He 

H 1 u f 

found the exponent m = 0.385 for all the weirs tested. 
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Thus the discharge equation for all sharp crested weirs 
was obtained by Villemonte as 

Q H, n 0.385 

7T = [ 1 - (—) ] (B.1.6) 

Q f H 1 

Further Villemonte recommend the following for practical 
application of the equation (B.1.6): 

(a) Measurement of upstream head should be done at a 
distance of at least 4H^ from the crest. 

(b) Measurement of downstream head H -2 should be done just 
beyond the surface turbulence caused by the nappe. 

(c) (i) For proportional weirs about five percent of 

P 

accuracy can be expected if — rr” > 1.0# 

H 1 

(ii) For rectangular weir three percent of accuracy can 

p 

be expected if — — >3.0 and 

H 1 

(iii) For triangular and parabolic weir two percent of 

p 

accuracy can be expected if — — > 3.0. 

H 1 

B.1.3 Comments on Villemonte's Work 

The following major comments can be offered on the 
Villemonte work. 

(1) Villemonte used the best fit line# drawn by eye 

H 2 Q 2 

judgement on the logarithmic graph 1 - (— ) n vs. rr- 
to obtain the value of m = 0.385 and K g = 1.0. The 
regression analysis of this data indicate slight variation 
in the value of m as given in Table (B.t.1) and hence the 

, suggested m = 0.385 is only an approximate value. 

II 
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Table B.1.1 

Value of m for K g = 1.0 for Villemonte Data 


SI. 

No. 

Type of weirs 

n 

m 

Percentage 
devi ation 
from 0. 385 

1. 

Linear proportional 
weir 

1.00 

0.368 

4.4% 

2. 

Rectangular weir 

1.50 

0.370 

3. 6% 

3. 

Parabolic weir 

2.00 

0.390 

1.8% 

4. 

Triangular weir 

2.50 

0.408 

5.9% 

5. 

Cusp parabolic weir 

3.32 

0.404 

4.9% 


(2) The range of the variables used by Villemonte 


are shown in Table (B. 1.2) 



Table B. 1* 2 

Range of Variables used by Vlllemonte's Study 
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From this it is seen that the range of — — is too 

H 1 

restrictive and studies of weir behaviour at lower value 

p 

of rr~ is desired. 

H 1 

(3) While from boundary condition K =1.0 for H,=0 

s 

it is possible that over a wide range of submergence ratio 
H 2 

rr- the equation (B.2.5) may have good correlation with 
K s / 1.0, i.e. Q = K [ 1 - . In such a 

XT S c li^ 

case determination of the best fit m and K g by regression 
analysis may give better accuracies. 


(4) Villemonte has used only one weir shape in the 
range n <1.5 viz. linear proportional weir. It is 
necessary to study the behavior of equation (B.-1.5)at more 
weir shapes with n <1.5. For example, in addition to 
sutro weir, the behavior of a quadratic weir could also be 
explored to generalise the use of equation B.5. 


B.1.4 The Present Study 

While Villemonte equation B.l. 6 is extensively quoted 
in many test books and hand books on fluid mechanics 
(Ref. 2,4, 5,7, 8), it is surprising that there is no other 
published literature on submerged flow over sharp crested 
weir. Also, some of the aspects of the flow such as the 
comments above, need to explored to validate the use of 

P 

Villemonte equation for a wide range of n-values and - 77 — 

rt l 

values. As such an experimental investigation of six 
different shapes of weirs were done in the hydraulics 
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laboratory of the Civil Engineering Department, Indian 

Institute of Technology, Kanpur. 

B.1.5 Objective of Present Study 

The main objectives of the investigation were: 

(1) to study the validity of Villemonte equation (B.1.6) 
for the range of . n from 0.5 to 3.5 and also over a 
wider range of submergence ratio. 

(2) to study the variation of the defined modularity 
limits as a function of n. 

(3) to note other aspects of submerged flow that may 
affects its use as an adequate flow measuring 


device 



CHAPTER B.2 


EXPERIMENTAL STUDY 


B.2.1 Introduction 

In the experimental study conducted at Hydraulic 
Laboratory# Indian Institute of Technology/ Kanpur. Six 
different type of weirs were used in free flow conditions 
as well as in submerged flow conditions. The weirs were 
selected in such a fashion that the exponent n in the 
discharge equation Q = K H^ varied from 0.5 to 3.5. 
For the present study the following weir shapes were 
used : 


Notional value of n 


(1) 

Quadratic weir 

0. 50 

(2) 

Linear proportional weir 

1.00 

(3) 

Rectangular weir 

1.50 

(4) 

Triangular 90° weir 

2.50 

(5) 

Parabolic weir 

2.00 

(6) 

Cusp parabolic weir 

3.34 


The photographic views of each of the above weirs 
are shown in Fig. B.2.1. 




••Jj 
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Fig . B. 2 . 1 . 3 


Fig. B.2.1.4 


Photograph of Weirs 




DIABOLIC WEIR,N*3-4 



B.2.2 Description of Experimental Set-up 


Channel : The experiments were conducted in a rectangular 
flume of length 4.5 m, width 0.45 m and depth 0.9 m. The 
dimensions and a line sketch of the flume are shown in 
Fig. B»2*,2. The water supply to the flume was from a cons- 
tant head overhead water tank. A valve in the supply 
line was used to control the flow. The water from the 
flume was diverted to sump through a calibrated measuring 
weir. 

Discharge Measurement: The water from the flume was passed 
through a linear ground level rectangular channel of width 
0.48m, which was fitted with a standard rectangular supp- 
ressed weir at the end of channel. This weir was well 
ventilated and carefully calibrated in this study by Volu- 
metric method. The calibration equation of this measure- 
ment weir was found to be 

Q = 1.048 h 1 * 5077 (B. 2.1) 

Q = discharge in m^sec,H = head measured at a 
specific stilling well in meters 

This calibration was valid for the range of head from 3 cm 
to 18 cm. The head on the weir was measured by a point 
gauge to an accuracy of 0.1 mm* 




Fig.B.2.2 Schematic sketch of the experimental layout 
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For a given discharge in the flume the head over the 
experimental weir as well as head over discharge measuring 
weir were taken with the help of point gauges. About 
25 to 40 sets to readings were taken in the practical range 
of discharge in the set-up. A total of 185 readings were 
taken in free-flow studies. 

After taking the readings in the free flow condition 
the weir was placed at a distance of 1.5 m from the supply 
tank inside the flume for submerged flow discharge studies. 
For a fixed discharge the desired submergence was produced 
with the help of the tail gate at the end of flume. The 
upstream head as well as the downstream head were measured 
at specified locations. In general a submergence range 
from five percent to 95 percent was produced and the corres- 
ponding head on the upstream and downstream positions were 
taken. A minimum of three sets of fixed discharge were 
used. The nature of the nappe in the downstream were also 
observed for different submergence ratios. About 30 to 45 
readings were taken for each type of weir# leading to a 
total of 225 reading for submerged condition. Basic data# 
viz.# H-l# H 2 and measured discharges are shown in Table 


B.2.1 



T ABLE B.E.1 

OBSERVED SUBMERGED FLOW DATA 


TABLE B . £ . 1 . 1 
QUADRATIC WEIR 

UPSTREAM HEAD DOWNSTREAM HEAD SUBMERSED FLOW DISCHARGE 

Ml (cm) HCicmJ Q5tlitr«*/s*c) 


34.41 000 

25.80000 

47.69910 

33 . 03000 

24 . 40000 

47 . 6771 0 

31 . 58000 

23 .30000 

47 . 67710 

£9.85000 

21 .80000 

47 . 67710 

ere . si ooo 

21 .20000 

47 . 6771 0 

£7.29000 

19.50000 

47.69910 

26.50000 

17.00000 

47 . 67710 

26 .2:0000 

1 5 . 50000 

47 . 6771 0 

2 4 . 57000 

1 T .70000 

47 . 6771 0 

22 . 50000 

1 5 . 76000 

47 . 67710 

c 0 . 95000 

13.25000 

47.69910 

zro . T^OOD 

13.20000 

47 . 6771 0 

1 5 . 3TDDO 

07 .46000 

47 . 6771 0 


OT . 1 0000 

4T . 6771 0 

33 . 93000 

21 . 66000 

57.86529 

27 . otooo 

1 5 . OOOOO 

37 . 56327 

25 . 13000 

1 T . 20000 

37 . 56327 

27 . 35000 

16.60000 

37.86327 

86.87000 

15.60000 

57.86529 

23 . 6T000 

13.60000 

ST . COSE? 

23.32000 

10.00000 

ST . BOSE? 

23 . 1 3000 

05 . 00000 

ST .B6SE7 

34.11 000 

£4 . 1 0000 

50.99486 

31 .01000 

22 .00000 

SO . ??4B6 

31 . 07000 

21 . TOOOO 

30 . 774B6 

27 . TCOOO 

20 . 50000 

SO . 774B6 

22T . 37000 

1 5 . 70000 

30 . 774B6 

87.38000 

18.80000 

50.99486 

25 .30000 

1 7.20000 

50 .99486 

24 . 75000 

1 6 .80000 

50 . 99486 

23 . 50000 

1 4 . 80000 

50.99486 

28 . 6 T 00 Q 

14.1 0000 

50 . 99486 

£ 1.60000 

12.36000 

50.99486 

21.12000 

1 1 .60000 

50 . 99486 

1 7 . 50000 

08 . 60000 

50 . 99486 

18.87000 

06 .70000 

50 . 99486 


112168 
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UPSTREAM HEAD 
Hi < cm ? 


4c . 45000 
CC . 05000 
C4 .55000 
53.5^000 
Cl . 37000 

El . 18000 

CO . 050CO 

17.70000 
Cl .04000 
£5 . 77000 

31 . 18000 

33.70000 
3T .75000 
40 . 75000 

48 . £7000 
44 . 44000 
4 T .0300 0 
34 . 00000 

30 . 1 5000 
31.1 COOO 

£4.70000 
EC . 554000 
CO . T5000 
17.13000 
1 5 . ©3000 

f 5 . 64000 
I 4 . 35000 
13 . 38000 
13 . 8C00O 

£5.70000 
E7 . 7EOOO 
30 . 06000 
EO . 33000 
El . 6EOOO 

££.86000 

£3.7£000 

18. 15000 
15 . £0000 
18.87000 


TABLE B. 8. 1 .8 

LINEAR PROPQFTI OMAL WE IK (5 *1 .01 

DOWNSTREAM HEAD 511 EMERGED FLOW DISCHARGE 

HE( cm i gbsc i itres/sec i 


35 .330 00 
CO . 530 OD 
1 3 . 530 00 
I . 1 ©0 00 
10.15000 

09 . 60000 
07 . 53000 
04 . 40000 

1 0 . 70000 
1 5 . £3000 

£4. 70000 
C5 . 50000 
33 . 50000 
43 . 30000 

40. 1OQO0 
33. 33000 
3 1 . 40000 
CO. 33000 
1 0 . 80000 

13. 20000 

£1 . £3000 

15.70000 

18.75000 
1 4. 7300 0 
11.80000 

05. 5800 0 
O7.E3OO0 

04 .00000 

03.75000 0 

£0.30000 
£3 .00000 
£4.80000 

13.20000 
14.85000 

16.60000 
17.60000 

08.90000 

1 1 .00000 

06 .90000 


54 . 18480 
34 . 1 8450 
54 . 1 3480 
34 . 1 8430 

34 . 1 5450 

54. 18480 
34 * 1 3430 
34 . 1 8460 

34 . 1 8480 
34 . 18430 

54.18480 

54 . 1 5450 
54 . 1 3430 

34 . 1 8480 

8£. 6E000 
5C. 3CO00 
6C . 8C0O0 
8£ . 6£0O0 
5C. ft CO D 0 
8C. ftCOOO 

34 . £5580 
34 . £5580 
34 . C5580 
34 . C5530 
34 .C5S30 

34. £5580 
34. £5580 
34. £5580 
34 .£5580 

44.70180 

44.70180 

44.70180 

44.70180 

44.70180 

44.70180 
44.70180 
44.70180 
44.70180 
50 . 55486 


TABLE B . 2 . 1 ,3 
RECTANGULAR WEIRCn«1 .5) 


UPSTREAM HEAD 
HI ( cm ) 


£ 1 . 85 0 0 0 
20.40000 
19.17000 

18.97000 

1 8 . 76000 

1 7 . 98000 
1 6 . 55000 
16.31000 
15.86000 

15.76000 

15.33000 
13. 09000 
13.00000 
13.41000 
1 3 . 39000 

18.61 000 
1 T . 39000 
1 T . COOOO 
16.00000 
1 D . 30000 

16 .20000 
1 3 . 93000 
1 3 . 47000 
1 3 . 44000 
1 4 . T! 000 

13.79000 

1 3 . 33000 
13 . OCOOO 
12.36000 
1 2 . 62000 

20. 19000 
1 9 . 20000 
1 T . 45000 
17.16000 
1 5 . 58000 

14 . 95000 
1 4 . 63000 
14.43000 
14.35000 
1 4 . 23000 

13.97000 
13.91000 
1 3 . ©4000 
13. ©0000 
1 3 , ©2000 


DOWNSTREAM HEAD SUBMERGED FLOW DISCHARGE 

H2 ( cm 1 QSC Iltr*£/S#c> 


18.90000 

1 7 .66000 ’ 

15.56000 

1 5 . 24000 

1 4 . 92000 

57.43930 
57 . 43930 

57.43930 

57 . 43930 

57 . 43930 

13.50000 

09.90000 

03.90000 

07 . 40000 

07 . 1 0000 

57.43930 
57 . 43930 

57 . 43930 

57 . 43930 
57.43930 

05.48000 

04.66000 

04 . 04000 

03.32000 

03 . 1 © 000 

57.43930 
5T .43930 

ST . 43930 

ST . 43930 

ST . 43930 

16.68000 

1 3 . 32000 

1 4 . T4000 

1 3 . T3000 

13.34000 

45.09300 

43 . 09300 

43 . 09300 

43 . 09300 

43 . 09300 

13.06000 

1 E . T6000 

1 1 . CEOOO 

1 1 . 56000 

1 0 . 44000 

45.09300 
45 . 09300 

45 . 09300 

45 . 09300 

45 . 09300 

07.64000 

05 . 55000 

0* . CEOOO 

04 . 05000 

OE .50000 

45.09300 
45 . 09300 

45 . 09300 

45 . 09300 
45 . 09300 

18.00000 

1 6 . 56000 

14.06000 

1 3 .66000 

05 .50000 

51 .93300 
51 . 93300 
51 .93300 
51 .93300 
51 . 93300 

07.80000 

07 . 04000 

06 . EEOOO 

05 . 5E000 

04 . 50000 

51 .93300 
51 . 93300 
51 .93300 
51 . 93300 
51 . 93300 

04.24000 

03 .T4000 

03 . E4000 

03 . EEOOO 

01 . T9 0 0 0 

51 .93300 
51 .93300 
51 .93300 
51 .93300 
51 . 93300 


TABLE 8.8. 1 .4 
PARABOLIC WEIR 
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UPSTREAM HEAD 
Ml tern? 


83 . 45000 
tl.07000 
ET . TOOOO 
C7 .©5000 
ET . 33000 

85 . 37000 
E4 . or c c <? 

E3 . T 4 U 0 0 
E3 .00000 
EE . ©3000 
EE . T3000 
EE . 05000 

30.67000 
3 t . 55000 
35.77000 
31.73000 

31 . 55000 
30 . 51000 
53.00000 
CT.TIOOO 

£7.68000 
ET . 5S000 
ET . 35000 
ET . 1 0000 
ET . 0E000 

£5 . 98000 
86 . 1 9000 
86 . 50000 
87 . 85000 
85 . 65000 

87.17000 

85.81000 
85 . 1 4000 
84 . 95000 

84.81000 


DOWNSTREAM HEAD SUBMERGED PLOW DISCHARGE 

HE C cm 3 G5(l ltrtl/ftc J 


18.80000 

13.00000 
E3 .00000 

56 .00000 
E3 . 30000 

17.90000 

13.50000 

13.30000 

10.00000 
OT . 00000 
04 . EOOOO 

01 . 90000 

£3.96000 

EG . 1 3000 
55 .EOOOO 

31 .40000 

86 . 00000 
E3 . 1 0000 
7 O . 63000 
7 4 . T3000 

14. 13000 
7 3 . 36000 
1 0 . T3000 

08 . 90000 

06 .90000 

15.83000 

1 6 . 40000 

17.43000 

1 9 . 80000 

83 . 43000 

19.53000 
75.1 0000 

1 1 .80000 
09 . 03000 
06.10000 


33 . 79669 
33 . T9669 
33 . T9669 
33 . T9 6 6 9 
33 . T9669 

33.79669 

33 . T9669 
33 . T9 6 69 
33 . T9669 
33 . T9669 
33 . T9669 
33 . T9669 

48.67400 

46 . 6T400 
48 . 6T400 
46 . 6T400 

48.67400 
48 . 5T400 
48 . 6T400 
48. 67400 

48.67400 
48.67400 
48.6T4O0 
48 . 6T400 
48 . 6T4O0 

40.35700 

40.35700 

40.35700 

40.35700 

40.35700 

40.35700 

40.35700 

40.35700 

40.35700 

40.35700 



TABLE B . £ . 1 .5 

TRIANGULAR 90-DEGREE WEIR 
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UPSTREAM HEAD 

Hi (cm) 

DOUNSTREAM HEAD 

HE ( cm ) 

SUBMERGED FLOU DISCHARGE 
QS ( l ltres/sec ) 

18.39000 

06 . 1 7000 

£0.04960 

18.51000 

08.17000 

20 . 04780 

18.61 000 

03 . 37000 

20 . 04780 

1 5 . 78000 

1 1 .53000 

20 . 04780 

1 7 . EOOOO 

12.81000 

20 . 04980 

1 9 . 47000 

13.95000 

£0.04980 

20.11000 

1 6 . 03000 

20 . 04780 

21 . 28000 

1 7 . 57000 

20 . 04780 

22 .50000 

20.55000 

20 . 04780 

23.58000 

21 . 85000 

20 . 04780 

19.77000 

15.27000 

£0.04980 

15.70000 

1 0 . 55000 

20 . 04780 

18.66000 

10.13000 

20 . 04780 

18.-10000 

OT .27000 

20 . 04780 

18.33000 

05 . 65000 

20 . 04780 

1 5 . »8000 

02.81000 

20 . 04780 

1 9 . 42000 

07.59000 

£3. 0£300 

1 5 . 3*3000 

06.07000 

23 . 02300 

1 3 . 28000 

04.31000 

23 . 02300 

1 7 . 32000 

04 . 73000 

23 . 02300 

17.36000 

05.71 000 

23 . 02300 

1 9 . 40000 

06.67000 

23.02300 

1 7 . 4T00Q 

08.25000 

23.02300 

1 7 . 5-300 0 

07.31000 

23 . 02300 

20 . 18000 

13.11 000 

23 . 02300 

20 . 42000 

1 4 . 15000 

23. 02300 

£ 0.97000 

16.09000 

£3. 0£300 

22 . 17000 

1 8 . 77000 

23. 02300 

23 . 22000 

20 . 87000 

23. 02300 

21.1 0000 

1 6 . 83000 

23. 02300 

21 .06000 

1 6 . 67000 

23.02300 
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CUSP PARABOLIC UEIR(n*3.4) 

UP STREAM HEAD DOWNSTREAM HEAD SUBMERGED FLOW DISCHARGE 

niicmi nctcm? 03UiTr**/ , seo 


eo . 660 00 

02 . 32000 

15.34790 

CD . 67000 

05 . 5C000 

15 . 34T90 

Cl .05000 

1 C . SCO 0 0 

15 . 34T50 

CC.360D0 

1 5 . 5C000 

15 . 34T9 0 

C0.6C0D0 

05 . 3C000 

15 . 34T3 0 

£0 . 78000 

09 . 1 £000 

15.34790 

co . ccooo 

03 . 3C000 

15 . 34T30 

C I . 03000 

1 C . 3C0 0 0 

15 . 3 J !T3 0 

Cl . 10000 

1 3 . PCOOO 

15 . 34T30 

Cl . 53000 

1 5 ,‘OCOOO 

15 . 34T90 

££. 00000 

1 6 .52000 

15.34790 

C3 . TOOOO 

CO . SCO 0 0 

15 . 34T9Q 

Cl . 5CU00 

1 5 . SCO 0 0 

15 . 34T9 0 

CO . 04000 

03 . SCOOO 

15 . 34T90 

Cl .41000 

1 s . ccooo 

15 . 34T30 

1 3 . SCQOO 

05 . TCOOO 

15 . 34T90 

19.67000 

09.62000 

12.51725 

Cl . 15000 

IT. 6C000 

1C.51TCS 

Cl . 50000 

1 5 . 4C0 00 

1C.51TCS 

ec . tcooo 

Cl . 3C000 

1C.31TC5 

• 13,31000 

1 C . CCOOO 

1C.51TCS 

19.89000 

1 £ . 02000 

12.51725 

CO . 38000 

15.1 COO0 

1C .51TC5 

Cl . COOOO 

1 6. TCOOO 

1C .51TCS 

1 3 . 5S000 

1 1 . TCOOO 

1C .S1TC5 
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Weirs i The weirs used for the study were prepared from 
3 mm thin aluminium plate by cutting and filling accord- 
ing to specified dimensions. The standard requirement of 
the sharp crested weirs were considered in the preparation 
of the plate. Pig. B. 2.3 shows the details of the various 
weir shapes used. 

The experiments were conducted on the six different 
type of weirs in both free flow conditions as well as in 
submerged flow modes. 

Head Measurement : The head over the experimental weir 

in the flume was measured with the help of a point gauge 
fitted with Vernier of least count 0.1 mm. The head on 
the downstream of this weir in submerged condition was 
measured with the help of six piezometer tapping placed 
in the downstream portion of the weir as shown in Fig. 
(B.2.2) and connected to a bank of manometers. The water 

manometer reading were taken to an accuracy of 0.2 mm. 
B.2.3 Experimental Procedure 

Experimental Procedure* For free flow and submerged 
flow conditions the weir under study was first placed 
at the end of flume. The flow from the weir had a free 
over flow ensuring full aeration and free flow mode for 
all discharges. For free flow studies/ the flow in the 
flume was controlled by a valve in the supply line. The 
head on the weir was measured at a distance of about four 
times the upstream head after the flow was stabilised. 









Fiq. B.2.3-5 90 Triangular Fig.B-2.33 Rectangular suppressec 

(n-2-5) weir ( n=1-5 ) 
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B.2.4 Study of nappe: The natnr-o ^ 4-v, . 

xne na ture of the nappe in the sub- 
merged flow was observed for all the weir shapes. It was 

found that the nappe could be classified into two types, 
viz. (i) plunging nappe, (ii ) and the surface nappe. 


In a plunging nappe the water surface plunges into 
the downstream flow creating a depression of water surface 
in the intermediate neighbourhood of the weir (Fig. (B.3.4 ) ) . 
In contrast, the water surface in the surface nappe is 
smooth and there is no appreciable depression of the water 
surface on the downstream side (Fig. (B.3.5)). 


The plunging nappe was observed at low submergences 
(J^/H.^ of about 0.4 0) for rectangular, parabolic, 90° tri- 
angular and the cusp parabolic weir. However, for proporT 
tional weir and quadratic weir, there was no plunging nappe 
for all the submergence values of H 2 /H^ when the water 
surface was above the rectangular base weir. For high 
submergence ratio the surface nappe was observed in all 
weir shapes. The transition from plunging to surface nappe 
occurred at approximately a value of (H-^/H.^^ as 9^ ven 
below: 


SI. NO. 

Type of Weirs 

n 

Transition 

'Wt 

1 . 

Quadratic weir 

0.50 

0.00 


2. 

Linear proportional 
weir 

1.00 

0.00 


3. 

Rectangular weir 

1.50 

0.35 


4. 

Parabolic weir 

2.00 

0.55 


5. 

Triangular 90°-weir 

2.50 

6.60 


6. 

Cusp parabolic weir 

3.40 

0.65 




63 


In a plunging nappe, the oscillation of water surface 
in form of crest and trough was observed. The oscillating 
waves were found stabilised in the downstream at about four 
to five times of Head H^. However, in suppressed rectan- 
gular weir, n = 1.5 the oscillating water surface was stab- 
lised only at larger distance of about 8 H^ • This is pro- 
bably due to absence of lateral expansion of the downstream 
flow in the weir. 



CHAPTER B.3 


RESULTS AND DISCUSSION S 


B.3.1 Analysis of Experimental Data: 

The free flow discharge equation for different 
type of weirs were found by regression analysis between 
upstream head above the crest level of weir and the 
measured discharge Q f in the equation. 

(1) Q f = K (B.3.1) 

for the rectangular, parabolic, 90°-triangular and cusp 
parabolic weirs. 

(2) Q f = K (H t - (B.3. 2) 

for the quadratic weir. 

(3) Q f = K (H x - |) n (B.3. 3) 

for the linear proportional weir# 

where a = depth of the rectangular base weir in meter, 

H x = head above the crest of weir on the upstream in meter 
and Q f = measured discharge in m 3 /sec. 

The actual value o^ K and n obtained by the regression 
analysis of the free flow data are shown in Table (B.3.1). 



Free Flow Discharge Equation of Weirs by Calibration 
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In submerged flow of a weir with head and H 2 on the 
upstream and downstream respectively, the hypothetical free 
flow discharge to any upstream head H x was calculated as 
follows. The equations used were 


( 1 ) 


Q f = K . hJ 


(B.3.1 ) 


for rectangular, parabolic, 90 -triangular and cusp para- 
bolic weirs 


( 2 ) 


Q f = K (H x 


2a. n 
3 1 


(B. 3.2) 


for quadratic weir and 


(3) 


Q f = K( Hl - §) n 


(B.3.3) 


for linear proportional weir. The proper values of K and 

n were taken from Table (B.3.1). The submerged free flow 

discharge Q g for each set of upstream head and down 

stream head K 0 was known from discharge measurement. 

Q s 

Fig. (B.3.1 ) show the variation of against for 

all the weir shapes. It is seen that for each weir shape 
the submerged flow discharge decreases, i.e. Q s /Q f decre- 
ases with an increase of Detail explanation of 

this is given in the next section (B.3»4 .) 

Q 


Defining Q 


, H = 1 - (H 2 /H 1 ) n the dis- 


charge ratio as in the form of generalised Villemonte 
equati on 

H, 


. r _ , 2.n-im 

= K s t 1 - ( h7> 1 


(B.1.5) 


can be written as 










.B.3-1.6 Variation of Qs/Qf with H 2 /H 1 for cusp parabolic weir 
(n-2-0) 


t o 


Q r = K S (B.3.4) 

B.3.2 Range of Variables in the Present Study 

The range of variable used in the persent study 
for the submerged flow are as shown in Table B.3.2. 

It may be mentioned that compared to Villemonte 

work (Table (B. 2. 2)) , the following points also have been 
incorporated in the present work. 

p 

(i) The variable ^ cover the lower range 1.0 not 
studied by Villemonte. 

(ii) The quadratic weir (n = 0.5) not studied by 
Villemonte extends the range of weirs with n <1.5. 

B.3.3 Calculation of the Value of K g and m 

Using the equation Q r = K g the value of m 

for K = 1.0 were calculated for each weir shape in the 
present weir study; for the whole range of submergences, 
by regression analysis. The value obtained for different 
weirs are shown in Table (B.3.3). 

It is seen that the value of m are different from 
the Villemonte of 0.385. For an observed head set of 
and H 2 the submerged flow discharge was calculated from 
the equation 

^ = K s [l - CB.1.5) 

Using the appropriate K g and m as in Table B.3.4 and the 
error prediction was found by comparing the estimated 
discharge Q g with the observed discharge. 


Variables of Present Study 
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It is seen that the Rms percentage error against 
observed discharge for quadratic and linear proportional 
weirs are 9.106 and 9.103 respectively. However for 
other weirs n 1.5 the Rms percentage errors are in the 
range 1.976 to 3.814. 

Next the Villemonte data along with the data of 
present study were analysed as one group, with K = 1. 0 
in the equation (B.3.4) and the value of m were calculated 
by regression analysis for each weir. (Table B.3.4) shows 
the result. 


Table B.3.4 

The value of m for all Range of Submergence for K = 1 

3 

For the Data (Present Study + Villemonte) 


SI. No. 

Type of weirs 

n 

m 

Coefficient 
of correla- 
tion 

1 . 

Quadratic weir 

0.5 

0.1534 

0.925 

2. 

Linear proportional 
weir 

1.0 

0.2955 

0.940 

3. 

Rectangular weir 

1.5 

0.3719 

0.991 

4. 

Parabolic weir 

2.0 

0.3660 

0.983 

5. 

Triangular 90°-weir 

2.5 

0.3940 

0.994 

6. 

Cusp parabolic weir 

3.4 

0.4029 

0.989 


6 
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By comparing the Table B.3.4 with the Table B.3.3 
it is seen that the value of exponent m slightly increased* 
The observed discharge was compared with submerged flow 
discharge calculated from equation (B.1.5) using appro- 
priate value of K g and m from the Table B.3.4 and the 
error prediction was found out. 

The Rms percentage error obtained against observed 
discharge for quadratic and sutro weirs are 8.82 and 9.68 
percent respectively. In the case of n > 1.5 the Rms 
percentage errors found within 3.75. Due to the higher 
value of errors against observed discharge for quadratic 
and linear proportional weir# it was thought that the sub- 
merged flow discharge should be divided into two zones of 
submergence ratio and the value of m and K should be 
obtained for better result. 

By trial and error different range of submergence 
ratios were selected to split up the data into two ranges 
of submergence. The value of K and m were obtained by 

o 

regression analysis in that range. For particular range 
of submergence error were calculated. 

The optimal two ranges of submergence were selected 
as: first zone 05 to 50 percent of submergence and second 

zone 4 0 to 95 percent submergence zone. The combined data 
of the present study and Villemonte data were considered in 
these two ranges. For submergence in second zone (40-95) 
percent the value of K s and m were calculated by regression 
analysis shown in Table B.3.5, 



Date (Present Study + Villemonte) 
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Since K s *= f (n) and m — f (n) had no good convenient corre- 
lation# a small change was made in K and m values to obtain 

•*> 

convenient correlation with n for easy practical use. Thus 
the selected value of K g and m ares 

Kg = 1.46 - 0.46 (n - 0.5) (B.3.5) 

for 0.5<n <1.5 
K = 1.0 n > 1 • 5 

S 

m = 0.37 

For 0.5<n<3.5 

In the Fig. B.3.2 the graphical representation of variation 

of value of K and m with n are shown. It is noted that 
s 

this result is not the optimal correlation equation but 
the results are convenient and adequate for practical use. 
The error in the prediction of discharge was found by the 
equation 




H 0 _ 

(_l) n 

H 1 



(B.l. 5) 


against observed value is shown in Table B.3.6. The value 

of K and m were selected from equation (B.3.5). The 
s 

minimum Rms percentage error was obtained for cusp parabolic 
weir (i.e. 2 percent) and the maximum Rms percentage error 
was for the quadratic weir (i.e. 5 percent). In Table B.3.8 
the Rms percentage error for each weir is shown. 
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TABLE 8 . 3 . f 

:RROR IN DISCHARGE E ST I HAT I ON IN ZONE II SUBMERGENCE RAT I 0 ( 4 0- 95 ) PER CENT 

n m 

*T EQUATION Q»/Qr-K»( 1-«HE/H1 ) ) 

Ki» 1 . 46- . 46 < n-0 . 5 > for 0.5< n >1.5 
K-s= 1 . 0 for n > 1 . 5 
in” 0 . 37 


TABLE B . 3 . 6 . 1 

QUADRATIC UEIR m* . 37 Ks*1.46 


( Q 5 ) 

( Qu ) 

(Qe) 

(Hr ) 

(P*) 

30.617 

33 . 635 

47 . 699 

. 73781 

-6.1178 

49 . 49 S 

32.906 

47 .699 

. 73032 

-3 . 771 0 

48 . 329 

32 . 1 24 

47 . 699 

. 73331 

-1 .3215 

47 . 948 

3 1 .906 

47 . 699 

. 71455 

- .5223 

46.916 

31 .214 

47 . 699 

.71 698 

1 . 6426 

47 0S0 

3 1 317 

47 . 699 

.70992 

1 .3603 

45 707 

30 . 424 

47 . 699 

. 70885 

4 . 1 763 

44 . 435 

29.611 

47 . 699 

. 68821 

6 .8431 

44.309 

29.943 

47 . 699 

. 63246 

6 . 0587 

44 . 295 

29 . 395 

47 , 699 

. 63645 

7 . 1 356 

46 . 023 

30.906 

47 . 699 

. 51497 

3.5107 

5 0 . £4 5 

40 . 247 

57 . 865 

. 63837 

-4 . 1 120 

56.161 

37 . 552 

57 . 865 

. 61920 

2.9451 

55 . 55 8 

37 . 169 

57 . 865 

.61101 

3.9707 

>6.174 

37.600 

57.865 

.59541 

2.9223 

>5 . 762 

37 . 348 

57 . 865 

.58057 

3.6347 

>6 963 

38 . 229 

57 . 865 

. 52980 

1 .5597 

55 . 435 

36 .903 

50 .995 

.70654 

-8.7069 

>3 78 0 

33 . 823 

50.995 

. 69598 

-5.4606 

53 . t 08 

33 373 

50 . 995 

. 69797 

-4.1447 

>1 .SS9 

34 . 337 

50 . 995 

. 69987 

-1 . 1 062 

(9.700 

33.117 

50 . 995 

. .69003 

2.5388 

19 .912 

33 . £64 

50 . 995 

.68663 

£ . 1235 

(7.919 

3 ( . 947 

50 . 995 

. 67984 

6 . 0321 

(7.315 

3 i . 546 

50.995 

. 67879 

7.2162 

(8.541 

32 441 

50 . 995 

. 62979 

4.8119 

(7 .825 

31 . 974 

50 . 995 

.62197 

6.2163 

[8.815 

32 7 06 

50 .995 

.57222 

4.2747 

[9.183 

32 . 984 

50 . 995 

.54924 

3 .5521 


DISCHARGE IN Litres/sec 

=es t i max ed submerged flow discharge for Ks=1.46 m*0.37 
=estimated submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

= e * per im e n t a 1 submerged flow discharge 
=submergence ratio 

^percentage error in submerged flow discharge estimation 
between Gts and Ge 
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TABLE B . 3 . 6 . E 
SUTRO WEIR M-.37 Ks«1.23 


(Q*s) 

( Qw ) 

(Qe) 

(HD 

(Pel 

S9 . £36 

46 . 484 

54 . 185 

.90209 

“9 . 3£1 6 

S5 . 84 3 

44 466 

54.185 

. 741 81 

-3.0605 

34.896 

43 . 9£4 

54.185 

. 64402 

-1.3119 

S4 . 754 

43 . 886 

54.185 

.60153 

-1 .0504 

S4 . 468 

43.803 

54 .185 

.50370 

- .5£35 

S4 . 4S£ 

43 . 434 

54.185 

. 70142 

- . -4940 

S7 . 3£0 

43 . 493 

54.185 

. 79217 

-5.7857 

36.133 

44 . 36£ 

54.185 

. 84366 

-3.5997 

37.980 

43.668 

54.185 

. 87460 

-7 . 0034 

6S.439 

65.227 

82.620 

.83074 

.£191 

81 . 507 

64 . 676 

62.620 

. 79500 

1 .347E 

81 . 670 

64 983 

82.620 

. 75426 

1 . 1493 

80 . 9E4 

64 . 853 

82 . 620 

. 60491 

2.0527 

37 . 397 

£9.368 

34.256 

.83851 

-9.1689 

35 . 945 

E8 . 474 

34.256 

.81 687 

-4.931 6 

35.567 

£8 . E57 

34 . 256 

.77801 

-3 . 8284 

35 . 07£ 

£8 . 003 

34.256 

. 68925 

"£.3841 

34.315 

£7 . 464 

34 .256 

. 6381 1 

-1734 

35.179 

£8 . E9 1 

34 . 256 

. 5 0 383 

-E.6965 

47.230 

37.491 

44.702 

.78988 

-5.6562 

47.210 

37 . 356 

44.702 

. 82973 

-5.6098 

44 . 805 

35.846 

44.702 

. 64929 

-.2406 

45 . 809 

36 . 582 

44 . 702 

. 68686 

-2.4758 

46.180 

36 804 

44 . 702 

.72616 

-3.3061 

46 . 930 

37 . 368 

44 . 702 

.74199 

-4 . 9844 

45.419 

36.442 

44 . 702 

.57292 

-1 .6042 


DISCHARGE IN Litres/sec 

Q5*est i mat td submerged flow discharge for Ks*l .23 m*o.37 
Qw=estimated submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

Qe*experimentdl submerged flow discharge 
Hr*submergence ratio 

Pe=percent age error in submerged flow discharge estimation 
between Qs and Qe 
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TABLE B . 3 . 6 . 3 

RECTANGULAR WEIR M* . 37 Ks-1.00 


c as ) 

C QW ) 

( Qe ) 

(Hr ) 

(Pe J 

55 . 356 

53 . 643 

57 . 439 

. 89245 

3. 62734 

55 . 676 

54.501 

57 .439 

. 86552 

2. 72177 

57 . 461 

56.316 

57.439 

. 81740 

-. 0371 1 

57 . 443 

56 . 335 

57.439 

. 80940 

-. 00575 

57 . 308 

56 . 235 

57.439 

.80165 

. 22828 

57 . 63£ 

56 . 708 

57 . 439 

. 75888 

-.33520 

59 . 868 

59 . 288 

57 . 439 

. 61224 

-4. 2281 9 

60 . 937 

60.443 

57 . 439 

.56180 

-6. 08968 

44 . 484 

43.225 

45 . 093 

. 89953 

1 . 35 01 4 

45 . 827 

44 . 688 

45 . 093 

.87182 

-1 . 6273 1 

44 . 4S0 

43 . 393 

45 . 093 

.66180 

1 .42559 

44 . 846 

43 . 887 

45 . 093 

.83605 

.5476 1 

44.700 

43 . 785 

45 . 093 

.82485 

.87124 

45.324 

44 . 437 

45 . 093 

.81310 

-.51 163 

45 . 032 

44.175 

45 . 093 

.80579 

.13568 

45 . 362 

44 .597 

45 . 093 

.77287 

- .5960 1 

44.961 

44.193 

45.093 

.77681 

.29169 

45.235 

44 . 599 

45 . 093 

. 721 1 0 

-.31549 

47 . 402 

47.002 

45 . 093 

.57262 

-5 . 12030 

50 . 975 

49.567 

51 . 933 

.89466 

1 .84472 

49 . 355 

45. 072 

51 . 933 

.88162 

4 . 96333 

50.589 

49 . 603 

51 . 933 

.81219 

2.58839 

50.206 

49 . 261 

51 . 933 

.80293 

3.32526 

53 . 0S0 

52 .507 

51 . 933 

. 64895 

-2.22803 

54.514 

54.104 

51 . 933 

.54019 

-4.96936 

54.152 

53 . 795 

51 . 933 

.50164 

-4.27247 


DISCHARGE IN Litres/sec 

Q s * e s t imated submerged flow discharge for Ks k 1.00 m * 0 . 37 
Qw^esti mated submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

Ge^ e x per t m ent a 1 submerged flow discharge 
Hr=submergence ratio 

pp^percentaqe error in submerged flow discharge estimation 
between Qs and Qe 



TABLE B . 3 . 6 . 4 


CQ5 ) 

PARABOLIC 
CQw ) 

WEIR M- 

(Qe ) 

. 37 Ks* 1 . 00 

(Hr ) 

(Pe) 

33 003 

32 839 

33.797 

. 54584 

2.34833 

32.360 

32. 121 

33.797 

. 65587 

4 . 19075 

32 612 

31 981 

33.797 

.85921 

3.50420 

34.312 

33 . 502 

33.797 

.89692 

-1 .52570 

32 . 624 

32 . 009 

33.797 

. 85361 

3.47016 

34.457 

34.118 

33.797 

.70556 

-1 . 95525 

32.857 

32 . 649 

33.797 

. 63202 

2.69026 

33.067 

32.071 

33.797 

.58551 

2. 16044 

46.674 

4b . 041 

48 . 674 

.78122 

4. 10982 

47.194 

46 . 431 

48 . 674 

.821 18 

3 . 04022 

47.461 

46 .561 

48 . 674 

.85480 

2.49204 

47.807 

46 . 079 

48 . 674 

.89894 

3 . 01458 

46.111 

43 . 336 

48 . 674 

. 82409 

5 . 26475 

46.438 

43 . 848 

48 . 674 

. 76389 

4 . 60526 

46, 466 

48.112 

48.674 

. 65701 

4.4951 1 

47.183 

46 . 903 

48.674 

. 53158 

3. 1 8736 

47. 286 

47 . 055 

48.674 

.51159 

2.85254 

39 . 59S 

39 . 658 

40 . 357 

. 58622 

1 . 14425 

39.464 

39.185 

40.357 

. 626 1 9 

2.21257 

39 . 449 

39 . 131 

40 . 357 

. 65774 

2. 25074 

39. 1 11 

38.690 

40 . 357 

. 72661 

3. 08838 

38.146 

37 . 539 

40.357 

.81780 

5.47820 

39. 207 

35 . 799 

40 . 357 

.71881 

2.84985 

39 . 382 

39 . 149 

40 . 357 

. 58504 

2.41598 


DISCHARGE IN Litres/sec 

Qs^eit i mated submerged flow discharge for Ks*1 .00 m* 0 . 37 
Qw=e 5 t i ma t e d submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

Ge*exper iroent al submerged flow discharge 
Hr*submergence ratio 

Pe -per ce nt age error in submerged flow discharge estimation 
between Qs and Qe 
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TABLE B . 3 . 6 . 5 
TRIANGULAR 90-DEGREE WEIR 


i Qs ) 

< Gw) 

< Ge? 

(Hr) 

(Pe) 

1 9 . 92 0 

1 9 . 866 

20 . 050 

. 50349 

. 64898 

1 9 . 966 

19.873 

£0 . 050 

. 60748 

.41770 

19.771 

19.646 

20 . 050 

.66719 

1 .39212 

J 9 . 659 

1 9 .502 

£0 . 050 

.71649 

1 .94724 

1 9 . 432 

1 9 . 200 

20 . 050 

. 79712 

3.08251 

19.799 

19.344 

20 . 050 

. 91333 

1 .25239 

1 5 . £26 

19. 025 

20 . 050 

. 77238 

4. 1 0052 

1 9.636 

19.562 

20 . 050 

.56578 

2.06228 

1 9 . 737 

1 9 . 672 

20 . 050 

. 54287 

1 .55762 

22.838 

22.706 

23. 023 

.64965 

. 80505 

22 . 776 

£2.613 

23 . 023 

. 69295 

1 . 07410 

£2. 607 

£2.375 

23. 023 

. 76729 

1 .80778 

22.581 

22.216 

23. 023 

.85566 

1 . 92019 

22 . 680 

£2 .207 

23 . 023 

.89879 

1 . 48872 

22.040 

£1 . 777 

23 . 023 

.79763 

4.26877 

22 . 097 

£1 . 840 

23 . 023 

.79250 

4. 02170 



TABLE B . 3 

.6.6 



CUSP PARABOLIC WEIR 

M- . 37 Ks* 1 

. 00 

( Q3 ) 

( Gw? 

(Qe? 

(Hr) 

(Pe) 

1 5 . 34£ 

15.300 

1 5 . 348 

. 59365 

. 03729 

15.125 

1 4 . 933 

1 5 . 348 

.83274 

1 .45217 

15.243 

IS .£03 

15 . 348 

.58583 

. 68269 

1 5 . 204 

15.135 

1 5 . 348 

.61706 

. 94022 

1 S . 378 

15.277 

15.348 

.73859 

- . 19449 

IS. 9S1 

IS. 840 

15.348 

.75091 

-3.93270 

14.995 

14.883 

15 . 348 

.761 99 

£.29808 

1 5 . 067 

14.983 

1 5 . 348 

.71088 

1 . 83345 

1 2 . 520 

12.377 

12.517 

.83270 

-. 02055 

1 2 . 626 

18 . 459 

12.517 

.85674 

-.86881 

1 2 . 723 

1 8 . 443 

12.517 

. 92776 

-1 . 64143 

12.480 

18 . 441 

12.517 

.61376 

.29866 

1 £ . 492 

18 . 455 

12.517 

.60432 

. 20391 

13 306 

1 3 . £03 

12.517 

.76835 

-6 . 30028 

1 2 . 432 

18.448 

12.517 

.59043 

.28518 


DISCHARGE IN Litres/sec 

Qc = est mated submerged flow discharge for Ks=1 .00 m= 0 . 37 
Gw=estimated submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

Qe=experimental submerged flow discharge 
Hr=sub mergence ratio 

p e =p»rce nt age error in submerged flow discharge estimation 
between Qs and Ge 
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Fig. B. 3-2 Variation of K s and m with n for submergence 
ratio 40 to 95 percent 
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For zone I with submergence five to fifty percent 
the value of K s and m were calculated for the best corre- 
lation in the equation (B.3.4) by regression analysis and 
values of K g and m obtained are shown in the Table B.3.7. 


Table B.3.7 

For Data ( Present Study + Villemonte) 

The value of K g and m for submergence (5-50 ) percent 


SI. No. 

Type of weirs 

n 

K s 

m 

Coefficient 
of corre- 
lation 

1. 

Quadratic weir 

0. 5 

1.073 

0.108 

0.850 

2. 

Linear proportional 
weir 

1.0 

1.096 

0.211 

0.800 

3. 

Rectangular weir 

1.5 

0.9907 

0.354 

0.970 

4. 

Parabolic weir 

2.0 

0.9844 

0.286 

0.970 

5. 

Triangular 90° 
weir 

2.5 

0.9889 

0.294 

0.980 

6. 

Cusp parabolic weir 

3.4 

1.004 

0.365 

0.988 


Since in the range of submergence ratio five to 
fifty percent the value of K = f(n) # and m = f (n) had 

no good correlation. So optimal value of K g and m as 
shown in Table B.3.7 were not used and for practical use 
the value of K 0 and m were taken as 
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K g = 1.08 for n < 1.0 

K s = 1.00 for n ^ 1.5 

m = 0.13 + 0.24 (n - 0.5) for n < 1..5 

(B.3.6) 


m * 0.37 for n > 1.5 

The graphical representation of the variation of the 
values of Kg and n are shown in Fig. (B.3.3). The Table 
B.3. 9 show the relationship between submergence ratio and 
percentage error in the submerged flow discharge estimation 
against observed discharge. 


The flow discharge estimation was done using 
equation 


~ = K s [l - (H 2 /H 1 ) n ] m (B.1.5) 

by using appropriate optimal value of K g and m from 
Eq. (B. 3.7)o 

It is seen that the maximum Rms error is three per- 
cent for linear proportional weir and minimum Rms error is 
for the cusp parabolic weir with a value of 0.4 percent. In 
the Table (B. 3.8) the Rms percentage error for each weir 
in the submergence ratio zone I of five to fifty percent are 
shown along with zone II for fifty to ninety five percent 
submergence ratio. 
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TABLE B.3.8 

ERROR IN DISCHARGE ESTIMATION IN ZONE I SUBMERGENCE RATI 0 < 5-50 IPERCEN 

n m 

8Y EQUATION Qs/Qf «Ks < 1 - ( H2/H 1 ) ) 

Ks« ! . OS for 0 . 5< n >1.0 
Ks* 1 . 0 for n> 1 . 5 
m* 0 . 1 3+ . £4 ( n- . 5 ) for n<1.5 
m*0 . 37 for n> 1 . 5 


TABLE B . 3 . 9 . 1 

QUADRATIC UEIR »*.13 K*«1.0B 


CQs ) 

(Qw ) 

(Qe ) 

(Hr ) 

(Pe ) 

46 . 7 32 

30.906 

47 . 699 

.51497 

£.02656 

AS . 449 

. 32 . 75£ 

47 . 699 

.40711 

2.62182 

37 . S4S 

40 .287 

57 . 865 

.41982 

. 55359 

37 . 632 

42 . 023 

57 . 865 

. 34387 

. 40380 

49 . 922 

34 . 528 

50 . 995 

.44103 

£.1 0450 

49.818 

33 . 940 

50 . 995 

. 36566 

£.30819 


TABLE B . 3 . 9 . £ 

LINEAR PROPORTIONAL UEIR m*.E5 Ks-1.08 


(Qs ) 

(Qw ) 

(Qe > 

(HD 

(Pe) 

52.300 

44 . 236 

54.165 

. 47499 

3.47772 

52 115 

44.31! 

54.185 

. 45455 

3.62081 

52 . 700 

45.619 

54.165 

. 37563 

2 . 73955 

53 . 642 

47.913 

54.155 

.22111 

1 . 00097 

52.188 

43 . 803 

54.185 

. 50370 

3 . 6B461 

78.314 

68 . 089 

82 . 620 

. 35557 

5.21202 

78 . 763 

6T . 470 

©2 . 620 

. 4241 6 

4 . 66817 

33 . 708 

26.291 

34 . 256 

. 50353 

1 . 59953 

34 . 100 

30.014 

34 . 256 

. 29595 

. 4558 0 

33 . 982 

29 . 268 

34 . 256 

. 39797 

. 80059 

43 . 674 

36 . 790 

44 . 702 

. 49036 

2 . 29972 

43.290 

37 . 850 

44 . 702 

.33195 

3. 15894 


DISCHARGE IN Litres/sec 
Qs=esi lraaTtd submerged flow discharge 
Ks*1.03 m=0.13 for quadratic weir 
Ks = 1 .08 m=0.2S for linear proportional weir 

Qw=est imated submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

Qe=e xper imental submerged flow discharge 
Hr=submergence ratio 

Pe =pe rcent age error in submerged flow discharge estimation 

between Qs and Qe 
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TABLE B.3.9 

. 3 



RECTANGULAR WEIR H 

*.37 Ks* 1 

. 00 

( Gs ) 

C GW ) 

( Qe ) 

(Hr ) 

(Pt) 

38 . 748 

30 41 3 

57 . 439 

.46658 

—2 . £7876 

38.701 

S3 . 387 

57 . 439 

. 45051 

-2 . 19650 

S8 . 840 

38 . 630 

57 . 439 

.35747 

—2 . 43863 

38.571 

38 . 407 

57 . 439 

.30881 

-1 . 96966 

38.835 

S3 . 7E4 

57 . 439 

.26933 

-2.46508 

38 . 739 

33 . 539 

57 . 439 

. 37768 

-2 . 35049 

39 . 333 

39.070 

57 . 439 

. 39641 

-3 . 27867 

43 . 948 

43.713 

45 . 093 

. 4 4 0 45 

-1 . 89703 

45 . 748 

43 .574 

45 . 093 

. 370E0 

-1 . 45173 

45 .877 

43 . 743 

45 . 093 

.31726 

-1 . 73937 

46 . 033 

43 . 938 

45 . 093 

. 221 87 

-2 . 08439 

51 .588 

51 .277 

51 . 933 

. 481 20 

. 66466 

SE . EE7 

51 .970 

51 . 933 

. 42956 

- . 56S55 

52.636 

52.423 

51 . 933 

. 38467 

-1 . 35287 

SE . 878 

52.700 

51 . 933 

. 34434 

-1 . 81878 

SE .E43 

32.101 

51 . 933 

.30351 

* . 59721 

52.535 

32.419 

51 . 933 

.26887 

-1 . 16007 

52.714 

32 . 620 

51 . 933 

.23410 

-1 . 50362 

52 . 496 

SE .404 

51 . 933 

.23333 

-1 . 08503 

52 . 820 

SE.78E 

51 . 933 

. 131 42 

-1 . 70879 



TABLE B.3.9. 

4 



PARABOLIC 

WEIR M* . 37 

Ks* 1 . 00 


( Qs ) 

(Qw ) 

(Ge) 

(HD 

(Pel 

33.691 

33 . 465 

33 . 797 

. 41852 

.31276 

34.211 

34.105 

33 . 797 

. 28790 

-1 . 22672 

34.603 

34 . 367 

33.797 

. 1 6268 

-£ . 38580 

34 .738 

34 . 731 

33.797 

. 05888 

-2.78473 

47 . 554 

47 . 085 

48 . 674 

.50074 

2 .30013 

47.990 

47 . 574 

48 . 674 

. 47294 

1 . 40548 

48 . 767 

48 . 507 

48 . 674 

.37869 

- . 19120 

48 . 686 

o 

in 

CO 

48 . 674 

.31321 

- . 02379 

49.186 

49 . 0 24 

48 . 674 

.23846 

-.92884 

40.172 

39 . 885 

40 . 357 

.431 9S 

.45949 

40.618 

40 . 437 

40 . 35? 

.34620 

-.64703 

41 .179 

41.100 

40 .357 

.22718 

-2 . 03631 


DISCHARGE IN Litres/sec 

Qs=estimated submerged flow discharge for Ks*1.00 m*.37 
Gw=estimated submerged flow discharge by Villemont* equation 
(Villemonxe value of Ks and m) 

Ge = e x per i me n t a 1 submerged flow discharge 
Hr=suomergence ratio 

p e = pe r c e nt a ge error in submerged flow discharge estimation 
between Gs and Qe 
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TABLE B.3.9.5 
TRIANGULAR 90-DECREE WEIR 


( Q3) 

(Qw) 

(Qe) 

(Hr) 

(Pe) 

£0.212 

20.171 

20 . 050 

.33551 

-.80966 

20 . 098 

20.012 

20 . 050 

.44138 

-.24031 

19.991 

19.866 

20 . 050 

.50349 

. 29307 

20.007 

19.943 

20.050 

.39620 

.21461 

20 . 1 22 

20 . 090 

20 . 050 

.30824 

-.35977 

20 . 250 

20.245 

20.050 

.15372 

-.99946 

23.093 

23.022 

23.023 

.39083 

-.30377 

23.157 

£3.118 

23.023 

.31489 

-.58259 

23.204 

23.189 

23. 023 

.22355 

-.78741 

£3.264 

23.242 

£3.023 

.25518 

-1 . 04847 

23.284 

23.251 

£3.023 

.29494 

-1 .13543 

23.240 

£3.190 

£3 . 0£3 

. 34381 

- ..9441 5 

£3 . 072 

£2.933 

£3.023 

.42373 

-.21136 

22.944 

£2.821 

23.023 

.47646 

.34444 



TABLE B.3.9 

. 6 



CUSP PARABOLIC 

UEIR 


( Qs ) 

(Qw) 

(Qe ) 

(Hr ) 

( Pe ) 

15.295 

15.295 

15.348 

. 1 1229 

.34142 

15.271 

15.261 

15.348 

.31513 

. 49844 

15.298 

15.293 

15.348 

.25725 

.32517 

15.274 

15.241 

15.348 

. 43838 

. 480 1 1 

1 5 .264 

15.220 

15.348 

. 47646 

. 54735 

12.479 

12.467 

12.517 

.34426 

.30929 

12.527 

12.458 

12.517 

.48907 

-.07885 


DISCHARGE IN Litres/sec 

Qs*tst tmat»d submerged flow discharge for Ks*1.00 m*.37 
Qw=estimated submerged flow discharge by Villemonte equation 
(Villemonte value of Ks and m) 

Qe=*experiment al submerged flow discharge 
Hr=3ubmergence ratio 

p e =p ? rcent age error in submerged flow discharge estimation 
between Qs and Qe 
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Table B.3.8 

Rms Percentage Error in Discharge Estimation 





Rms of 

percentage error 

SI. No. 

Type of Weirs 

n 

Range of 

submergence ratio 

JLi 

Zone I 

Zone II 




(5-50) percent (40-90) percent 

1. 

Quadratic weir 

0 o 5 

1. 9 

5.0 

2. 

Linear propor- 
tional weir 

1.0 

4-.1 

3.1 

3. 

Rectangul ar 
weir 

1.5 

1.8 

00 

• 

CM 

4. 

Parabolic weir 

K> 

• 

O 

1.5 

3.2 

5. 

Triangular 90° 

2.5 

0.7 

2.3 


weir 




6 0 

Cusp parabolic 
weir 

3.4 

0.4 

2.2 


From Table B.3.8 it is seen that all the weirs could 
be used in submerged condition with percentage error less 
than five percent. The different type of weirs tested can 
be arranged in the following order in accordance with their 
accuracy in discharge measurement in submerged modes 

(1) Cusp parabolic weir 

(2) Triangular weir 

(3) Parabolic weir 

(4) Rectangular weir 



(5) Linear proportional weir 

(6) Quadratic weir. 

It is interesting to note that the cusp parabolic and 
triangular weir can be used for discharge measurement in sub- 
merged condition upto seventy percent submergence ratio with 
an error within one percent. 

B.3.4 The Study o f Modularity Limit: 

The variation of the observed values Q with H 0 /BL 

r z i 

for various weir shapes is shown in Fig. (B.3.1). It is 
usual to define a modular limit of flow measuring device as 
the maximum ratio of tailwater head to the upstream head 
(both measured above the common datum, usual Iv the crest) 
at which flow is not affected. To give a qualitative value 
of the head ratio B^/B^ at which significant changes in 
the discharge Q f would take place. The modular limit M^ 
and M l0 are defined as below: 

M^ = Modular limit at five percent level: It is 
defined as the value of submergence ratio at which 

the real discharge in submerged mo^e deviates by five 
percent from the discharge calculated by the head discharge 
equation of free flow. 

M^q = Modular limit at ten percent: It is defined 
as the value of submergence ratio B^/B^ at which the real 
discharge in submerged mode deviates by ten percent from 
the discharge calculated by free mode discharge equation. 

The value of M 5 and M^ for the various weir shapes 
obtained by interpolation from the plot of Q r vs. 

(i.e. Fig. (B.3.1)) are shown in Table B.3.10: 



94 


Table B.3,10 

The Value of and M^q for Each Weir 


SI •No* 

Type of weir 

n 

M 5 

M 

10 

1 . 

Quadratic weir 

0.5 

0.40 

0o 5 5 

2 . 

Linear proportional 
weir 

1.0 

0.40 

0.55 

3 . 

Rectangular weir 

1.5 

0.25 

0.40 

4 . 

Parabolic weir 

2.0 

0.36 

0.50 

5 0 

Triangular 90° weir 

2.5 

0.44 

0.56 

6. 

Cusp parabolic 

3.4 

0.54 

0.66 


It is seen that (Table B.3.10) for the quadratic 

and sutro weirs the value of M 5 and M ^ 0 are constant at 
0.40 and 0.55 respectively. But for rectangular, parabolic, 
90°-triangular and cusp parabolic weirs (i.e. n > 1.5) the 
modular limit (M) varies as function of n. 

Considering M = f (n) = a(n-1.5)^ + C (B.3.8) 

By regression analysis the value of a and b was found out 
in equation (B.3.8 ) taking the initial value of M at n=1.5. 
So the value of M 5 and M 1Q can be expressed as 







M 5 - 0.185 ( n -l . 5) * 727 + 0.25 (B.3.9) 

For n < 1 . 5 

M 10 * 0.163 (n-i.5)' 715 + 0.40 (B.3.10) 

For n < 1*5 

The value of and M^ 0 can also be calculated by equation 

= K g [ 1 - (H 2 /ta 1 ) n ] m (B.1.5) 

by selecting proper value of K_ and m. However the equations 

(B.3.9) and (B.3.10) is believed to be useful in quick esti- 

mation of the possible submergence effect. 



CHAPTER 4 


CONCLUSIONS AND RECOMMENDATIONS 


4.1 Conclusions 

Based on an experimental study of the submerged 
flow in six kinds of weir shape, the following signifi- 
cant conclusions are obtained. 

In the generalised form of Villemonte equation 

C> H 0 „ m 

S _ _ r* *i / r N HI / m . 

oj = K s t 1 - { S7 ! 3 

the variation of the value of K,. and m were obtained for 
all the weir shapes. 

(1) For all the weir shapes tested, in general the value 

of K s and m, obtained from the regression analysis of the 
data, are different from Kg = 1*0 and ra = 0.385 suggested 

by Villemonte. 

(2) (a) If K s =1.0 is taken, then for weirs with n > 1.5 
(viz., rectangular, triangular, parabolic and cusp para- 
bolic) the corresponding value of m is very near 0.385. 

In other words, the Villemonte equation can be used in 
these weir shapes to predict the discharge within an 
error of about five percent. 

(b) However, if Kg = 1.0 is taken for weirs with n<1.5 
(i.e. linear proportional and quadratic weirs) the experi- 
mental data correlation gives m = 0.15, for quadratic 
weirs and m = 0.26 for linear (sutro) weir against the 

Villemont's recommended value of 0.385. As such, the use 





of Villemonte equation in these weirs may lead to large 
errors# upto as much as 25%. 

(3) For convenience of discharge prediction two zones 
of submergence as below are suggested! 

Zone I: Submergence Jf to 50 percent 
Zone II: Submergence 40 to 95 percent. 

For zone I # based on the experimental data# the following 
values of K„ and m are suggested in the generalised form 
of Villemonte equation (B.2.5) 

For submergence (iT-50) percent 

K = 1.08 for 0.5 <n <1.0 

s 

K s = 1.00 for n ^ 1.5 

and m = 0.13 + 0.24^n-0.5) (B.3.7) 

These values of K_ and m if used in the equation (B.2.5) 

s 

prediction of discharge with an accuracy of 97 percent 
for all weirs can be expected. 

For zone II# based on the experimental data# the following 

values of K_ and m are suggested for the equation (B.2.5) 

For submergence 40 to 95 percent 

K g = 1.0 for n > 1. 5 

K = 1.46-0.46 (n-0.5) for 0.5in^l.5 

s (B.3.5) 

and m = 0.37 for 0.5 < n < 3.5 



99 


If these values of K and m are used in equation 
(B.2.5) for submerged flow discharge estimation discharges 
with an accuracy of ninety five percent for all weirs can 
be expected* 

(4) The modular limit at five percent level (M g ) and 
ten percent level (M^q) are studied and their values are 
obtained as follows: 

Mg = 0.40 for 0.5 ^ n ^ 1.0 

M 5 = 0.185 (n-1.5) °* 727 +0.25 (for n > 1.5) 

(B.3. 9) 

and M-^q = 0.55 for 0.5ZnZl.0 

M 10 = 0.163 (n - 1.5)°* 715 + 0.40 for n^. 1.5 

(B.3.10) 

(5) The triangular and cusp parabolic weirs are the best 
for discharge measurement in submerged flow conditions and 
in submergences upto 70 per-cent# the discharge can be 
measured by these weirs with an accuraay of 99 percent. The 
order of listing of weirs in term of their accuracy in 
discharge measurement in submerged mode are as follows* 

Weirs Percentage Rms error 

(1) Cusp parabolic weir 1.8 

(2) Triangular weir 2.0 

(3) Parabolic weir 3.0 

(4) Rectangular weir 


3.0 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


?.°.A 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

5 

1 

1 

1 


1 

5 

1 

1 

1 

1 

1 

1 

1 

i 

1 


1 


ev AMOO KUMAR 


A PROGRAM FOR SELECTION OF DISCHARGE MEASURING DEVICES. 

THIS PROGRAM HELPS IN IDENTIFYING A SET OF POSSIBLE DISCHARGE 
MEASSURING DEVICES FOR FLOU MEASUREMENT IN A GIVEN SITUATIONS 
TWENTY SIX DEVICES ARE CONSIDERED. 


INTEGER SO 

OPEN I UN I T=3c' , FILE= ' SELECT . OUT ' ) 

WRITE I 6 , 786 ) 

FORMAT < 5 1 X , ' WELCOME FOR THE SELECTION OF'/IH , 

1 Ox , , 'RELEVANT STRUCTURES FOR'/IH , 

IOX, 'QISCHRGE MEASUREMENT WITH OR'/ 

1 Is, 'WITHOUT DISHARGE REGULATION'/) 

PAUSE 

WRITE I 6, SOS) 

FORMAT! 1 OX, 'LIST OF DISCHARGE MEASURING STRUCTURES CONSIDERED 
I OX , 'FOR THE SELECTION OF RELEVANT STRUCTURE OR SECTION'/ 

1 OX, 'LIST OF WEIRS'/ 

10X,'U> ROUND NOSE HORIZONTAL BROAD-CRESTED WEIR'/ 

10X, ' £ £ ) ROMIJIN MOVABLE MEASUR I NG/REGUL AT ING WEIR'/ 

1 OX, 'REF SOS,M.G, DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 OX, 'NETHEERLAND 1978'/ 

IOX, MS) TRIANGULAR BROAD-CRESTED WEIR'/ 

1 OX, ' £ 4 ) SROAD-CRESTED RECTANGULAR PROFILE WEIR'/ 

I OX, ' IS) FA I YUM WEIR '/ 

1 OX, 'REF BOS, M.G, DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 OX, 'NETHEERLAND 157S'/ 

10X,' (6) RECTANGULAR SHARP-CRESTED WEIR'/ 

1 OX , ' ! 7 ) V-NOTCH SHARP-CRESTED WEIR'/ 

1 OX , ' £ 8 ) CIRCULAR WEIR'/ 

1 OX , ' ( 9 ) PROPORTIONAL WEIR'/) 

PAUSE 

WRITE! 6, 889 ) 

FORMAT! J OX, '( t 0 )WEIR SILL WITH RECTANGULAR CONTROL SECTION'/ 

1 OX, ' £ 1 t JV-NOTCH WEIR SILL'/ 

fOX,'IJE)TR! ANGULAR PROFILE TWO-O I MENS 1 ON WEIR'/ 

10X, ' £ 1 3 ) TR I ANGULAR PROFILE FLAT-VEE WEIR'/ 

1 OX , ' I 1 4 ) 8UT CHER MOVABLE STANDING WAVE WEIR'/ 

10X, ' REF BOS M G DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 OX, 'NETHEERLAND 1578'/ 

1 OX, ' £ 15 IUES-STANDARD SPILWAY'/ 

I OX , ' REF 60S M G DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 OX, 'NETHEERLAND 1978'/ 

1 OX , ' ! 1 A ) CYL I NOR I CAL CRESTED WEIR'/ 

1 OX, 'REF BOS M G DISCHARGE MEASURING STRUCTURE ILRI'/ 

5 OX , 'NETHEERLAND 1978'/) 

PAUSE 

WRITE (6, 686) | 

FORMAT! 1 OX, 'LIST OF ORIFICE'/ 

1 0 X , ' £ 1 ) RECTANGULAR SHARP-EDGED ORIFICE'/ 

1 0 X , ' ! t' ) RADIAL OR TAINTER GATE'/ 

1 0 X , ' £ 3 ) CRUMP-DE GRUYTER ADJUSTABLE ORIFICE'/ 
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1 1 OX, 'REF 60S M G DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 t OX, 'NETHEERLAND 1978'/ 

1 1 OX, ' (4) METER GATE'/ 

1 10X,'(S) NEYRPIC MODULES'/ 

1 1 OX , ' REF BOS n G DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 1 OX, 'NETHEERLAND 1978'/ 

I 1 OX, 'LIST OF FLUMES'/ 

1 10X,'(1> LONG THROATED FLUMES'/ 

J I 0 X , '( £ ) MARSHAL FLUME'/ 

1 10X, '(3) H-FLUME'/ 

I ! OX, 'MI THROATLESS-FLUME '/ 

1 1 OX , ' REF EOS M G DISCHARGE MEASURING STRUCTURE ILRI'/ 

1 10X, 'NETHEERLANO 1978'/) 

PAUSE 

C ******9*** ************************* ******9*iit?r**#i'#*m*ii*#****9 

SiCt UR I TE ( 6 , 7 C 0 ) 

700 FORMAT ! 1 OX , ' 00 YOU WANT DISCHARGE REGULATION ?') 

WRITE! £, TOE > 

70E FORMAT! 1 OX , 'YES' , ' THEN') 

UR I TE (6,704 ) 

704 FORMAT! 1 OX, 'WRITE 1') 

UR I TE ( 6 , 708 ) 

706 FORMAT( I OX, ' IF NOT V 1 OX , ' WRI TE £'/) 

READ ( 6 , *>1PFS 
PAUSE 

WRITE! 6, 71 0) 

710 FORMAT ( 1 OX, ' DO YOU WANT TO PASS SEDIMENT OR TRASH ?'/ 

1 1 OX , ' YES ' , ' THEN'/IOX, 'WRITE 3'/ 

1 1 OX , ' I F NOT'/I OX, 'WRITE 4V) 

READ i 6 , * ) ISEO 
PAUSE 

UR I TEC 6 , 71 £ ) 

71 £ FORMAT ( 1 OX , ' 00 YOU WANT TO PASS TRASH ? ' / 1 OX , ' YES ' , ' THEN'/ 

1 10X, 'WRITE S ' / 1 0 X , 'IF DO NOT ' /I OX, ' WRITE 6'/) 

READ < 6 , * ) I TRS 
PAUSE 

UR I TE ( 6 , 90S I 

90 S FORMAT ( 5 OX , 'ENTER Q 1 MAX QIMIN Y 1 MAX Y 1 MIN ' ) 

UR I TEC 8,9 04) 

904 FORMAT! 10X, ' Q! MAX THE MAXIMUM DISCHARGE IN CUBIC METERS/SEC ' ) 

1 10X, ‘QIMIN THE MINIMUM DISCHARGE IN CUBIC METERS/SEC'/ 

1 1 OX , ' Y 1 MAX THE MAXIMUM DEPTH OF WATER IN METERS'/ 

1 10 X, ‘Y1MIN THE MINIMUM DEPTH OF WATER IN METERS'/) 

READ! 6,* )Q1 MAX, QIMIN, Y I MAX, Y1 MIN 
PAUSE 

URITECb , 906 ) 

966 FORMAT! 5 OX, 'ENTER P QHMX $ 6 ') 

UR I TEC 8 ,9 fifei 

905 FORMAT! 1 OX, 'P HEIGHT OF CREST IN METERS'/ 

1 1 0 x , ' DHMX MAXIMUM HEAD LOSS IN METERS'/ 

1 1 OX , ' B BREATH OF BOTTOM OF CONTROL SECTION IN METERS'/) 

READ C 8 , * )P , DHMX, B 
PAUSE 

IF! CPFS. EQ.£>WRITE(6, 776) 

1 F ( IPFS . EG . 1 JURITE C 6 , 77E ) 
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77* FORMAT! JOX, 'DISCHARGE MEASURING STRUCTRES') 

771 FORMAT! )0X, 'DISCHARGE REGULATING STRUCTURES') 

IF! ISED.EQ. 3)WRITE( 6, 774) 

IF USED. ED. -4) WRITE (6, 773) 

TFfTTRS.EG. 5) WRITE ( 6 , 77£) 

IF £ I TR3. ED. 6) WRITE !6, 771 ) 

774 FORMAT ( 1 OX , 'SEDIMENT OR TRASH TO BE PASSED') 

773 FORMAT!) OX, 'SEDIMENT OR TRASH NOT TO BE PASSED') 

77c FORMAT! 1 OX , 'TRASH TO BE PASSED') 

77) FORMATOOX, 'TRASH NOT TO BE PASSED') 

C ****4'*******#i*4#5t***44 , **'t**i*4*^*f^**.f^**** ***»***********< 

C IN PUT DATA REGISTRATION 

WRITE! 6,777 )YtMAX, Yt MIN, Qt MAX, G ? M I N , P , B , DHMX 
777 FORMAT! ) OX, '!) >Y) MAX =',£X,F)O.S, )X, 'METERS' 

1 / f OX , ' ! E ) Y1MIN «' , EX , F I 0 . 5 , IX, 'METERS' 

1 /) OX , ' ! 3 ) 0) MAX =' , EX, FI 0 .5, )X, 'CUBIC METER/SEC ' 

1 / 1 0 X , ' (43 GIMIN =' ,£X,FJ 0 .5, IX, 'CUBIC METER/SEC' 

1 / 1 OX , ' ! S ) P =' , £X,F) 0 . S, )X, 'METERS' 

1 /!0X,'<6) B = ' , £X , F 1 0 . 5 , IX, 'METERS' 

) / ) 0 X . ' ! 7 ) DHMX =' ,2X,F) 0 .S, IX, 'METERS'//) 

WRITE! 6 , n n 

)1) FORMAT! ) OX, 'ARE YOU SATISFIED WITH INPUT DATA'/ 

1 5 0 X , 'YES' , ' THEN'/T OX, 'WRITE !'/ 

1 1 0 X , 'IF NOT'/) OX, 'WRITE £'/) 

READ ( 6 , * ) SO 
lFISD.EQ.S) GO TO 222 

C. ***»*»*****»*.•***♦*•*$■.***.********************************** ***i 

H) MAX*Y) MAX-P 
H) M) N=Y1 M I N-P 
&AMA“Q) MAX/G) MIN 
write!*,*) ' GAMA = ' , GAMA 
QPN=Q1 MAX/S 
H£MAX*H 1 MAX-OHMX 
AML=H£MAX/H1 MAX 

IF ( KEMAX . LT . 0 . 000 I ) WRI TE ( 6 , 770 ) 

IF! H2MAX . GT .0.000) ) WR I TE ! 6 , 670 ) 

WRITE! *,£00) 

200 FORMAT! ) Ox, “DO YOU WANT TO GET RESULT IN A FILE NAMED ", 

I "SELECT . OUT "/ I Ox » " IF YES THEN WRITE 1 ”/1 Ox , "OTHERWISE £" } 
READ! *,*) IANS 
(F! IANS. EG. I ) GO TO £01 

GO TO £0£ 

£01 WRITE! 3£ , 904 ) 

WRITE! 32, SOS ) 

WR I TE ! 32 , 777 ) Y 1 MAX , Y 1 H I N , Q ! MAX , G 1 HI N , P , B , DHMX 
IF ! TPFS .EQ . 2) WRITE! 32, 776 ) 

IF! IPFS.EG. 1 ) WRITE! 3E , 775 J 
!F( ISED.EQ. 3)URITE(32,774) 

IF! ISED. EG. 4) WRITE! 3£, 773) 

IF! ITRS .EQ.S) WRITE! 32, 772) 

IF(T7RS.EQ.6)WRITE(3£,7?1 ) 

I F ! H2MAX . LT .0.0001 ) WR I TE ! 32 , 770 ) 

I F ( He MAX . GT . 0 . 000 I ) WR I TE ( 3£ , &70 ) 

770 FORMAT! ) OX, "FREE FLOW CONDITION"//) 

670 FORMAT! 1 0X , "SUBMERGED FLOW CONDITION”//) | 
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£0£ IT( jrPPS.EO. 1 >G0 TO 6 

1FCTPFS. £Q.£)GO TO . 5 

5 IFC ISED.ECt.3)G0 TO 1234 

IFC ISED.EG.4)G0 TO 1236 

1034 CALL SDPSCH1MAX,H£MAX,H1MIN,Q1MAX,Q1MIN,P,B, IANS) 

GO TO 92 

3 236 CALL SONOC HI MAX , H2MAX , HI M IN, Q1 MAX , Q 3 M IN , GAMA , I ANS ) 

GO TO 9£ 

6 i F C GAMA . GT . 35 . 0 > SO TO £S 
3 F £ GAMA . LE .35.0) SO TO £9 

£9 1FCITRS.EG.S) GO TO 1£35 

I F< ITRS . £0 . 6 ) GO TO 1£37 

1 £33 CALL TRPSC HI MAX, H£MAX, Gt MAX, Gt MI N,P,B, ISED, IANS) 

GO TO 92 

1 £37 CALL TRNO ( H 5 MAX , H2MAX , H 1 MI N , Q f MAX ,QfMIN,P,B,X SED , IANS } 

GO TO 9£ 

£0 I F £ GAMA . GT . 1 £0 )G0 TO 3 0 

IF (H1MAX.LT .£.0)60 TO 3£ 

J F ( I ANS .£0.1 >G0 TO 355 
GO TO 356 

355 UR I TE { 3£ , 59 ) 

356 WR I TE (6,59) 

59 FORMAT ( 1 OX.. “THE RELEVANT STRUCTURE IS"/ 

3 3 OX/ ” CYL ( NOE R I CAL CRESTED WEIR”/) 

GO TO 9£ 

3£ IFCHIMAX.LT. 1 . 0)GO TO 34 

IFC IANS .EG . 1 ) GO TO 357 
GO TO 350 

35? UEIT£C3£,61) 

GO TO 92 

358 WR I TE ( 6 , 6 3 ) 

61 FORMATC 1 OX, * , THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, " ( 3 ) CYLINDERICAL CRESTED WEIR”/ 

1 1 OX, “ (2! TRIANGULAR PROFILE TUO DIMENSIONAL UEIR*‘//> 

GO TO 92 

34 IFC IANS. EG. 1 >G0 TO 359 

GO TO 380 

359 WRITE (32, 83 3 

38 0 UR I TEC 6, £31 

61 FORMAT (t OX, “THE RELEVANT STRUCTURES ARE*'/ 

1 1GX,“C1> CYLINDERICAL CRESTED WEIR"/ 

3 3 OX , " ( 2 > TK I ANGULAR PROFILE TWO DIMENSIONAL WEIR"/ 

1 1 OX , “ C 3 ) BUT CHER MOVABLE UEIR“//> 

GO TO 92 

30 I F < H t MAX . LT . £ . 0 ) GO TO 36 

IFC IANS.EQ. 1 )G0 TO 363 
GO TO 362 

363 WRITE (32,65) 

3C £ UR I TEC £,651 

85 FORMAT C 3 OX , ’’THE RELEVANT STRUCTURE IS”/ 

1 1 OX, “CYLINDRICAL CRESTED UEIR WITH TUO OR MORE STRUCTURES"//) 

GO TO 92 

36 IFC IANS.EQ. 1 ) 60 TO 363 

GO TO 364 

363 WRITE (32, 67) 
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UR I TE t 6 , 6? ) 

FORMAT ( 1 OX, "THE RELEVANT STRUCTURES ARE”/ 

1 OX, *! 1 > CYLINDRICAL CRESTED WEIR**/ 

10X,"{E)TR[ ANGULAR PROFILE TWO OIHENSI0NAL WEIR", 
* WITH TWO OR MORE STRUCTURES " // ) 

CONTINUE 

STOP 

END 


SUBROUTINE TRNO ! HI MAX , HEMAX , HI M 3 N, Gt MAX, Q1 M I N, P , B , I3ED, IANS) 
PAUSE 

WRITEI6, III) 

FORMAT! 1 OX, 'ENTER FLX FLEXIBILITY'/) 

READ 16,*) FLX 

WRITE! 6, SSSIFLX 

FORMAT ! 1 OX, M S)FLX= ' , £X , FA . £) 

IFOANS.EQ. 1 ) GO TO £01 
GO TO £0£ 

WR I t c. ( 3c.' , 'zi'jj'zi )t-LX 
I F ! FLX . GE . I JCO TO 4E 
!F(H1MAX.LT.E. OGO TO 44 
WRITE! 6, 69) 

J F ! IANS . EQ . 1 ) GO TO £03 
GO TO 94 
WRITE! 32 , 69 > 

FORMAT! I OX, "THE RELEVANT STRUCTURE IS"/ 

1 OX, “MOVABLE CYLINDRICAL CRESTED WEIR 1 *//) 

GO TO 94 

IF ( HI MAX . LT . I . 0 ) GO TO 46 
WRITE (6, 7 I ) 

IF! IANS. EQ. 1 )GO TO 300 
GO TO 94 
WRITE! 3£, 71 > 

FORMAT (1 OX, "THE RELEVANT STRUCTURES ARE"/ 

1 OX, “ ! 1 ) MOVABLE ROUND NOSE HORIZONTAL BROAD CRESTED UEIR“/ 
10X, "iSIMOVASLE TRIANGULAR BROAO CRESTED WEIR”/ 

1 OX , “ ! 3 ) MOVABLE CYLINDRICAL CRESTED WEIR'V/) 

GO TO 94 

AML=H£MAX/H I MAX 

IF ( AML . LE . 0 . 33 ) GO TO 48 

WR ITE ! 6 , 71 ) 

3 F ! IANS .EQ. 1 )GO TO 301 
GO TO 94 
UR ITE! 32 , 73 ) 

FORMAT! I OX , "THE RELEVANT STRUCTURES ARE"/ 

1 0 X , " ! 1 ) MOVABLE ROUND NOSE HORIZONTAL BROAD CRESTED WEIR,”/ 

1 0 X , ” ( £ 3 MOVABLE TRIANGULAR BROAD CRESTED WEIR,”/ 

1 OX, “ ! 3 ) BUTCHER MOVABLE UEIR'V/1 OX, “ ! 4 ) DIVISORS"//) 

GO TO 94 
WR I T£ ! 6 , 7S ) 

IF! JANS .EQ .1)60 TO 302 
GO TO 94 
I Tc { 32 7G ) 

FORMAT ( 1 OX , "THE RELEVANT STRUCTURES ARE MOVABLE WEIRS’ 1 / 

1 OX, “ ! 1 )ROUND NOSE HORIZONTAL BROAD CRESTED WEIR,"/ 
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1 JOX, " ( £ ) HAM l JI N MOVABLE MEASURING OR REGULATING WEIR,”/ 
1 10X, ** £ 3 ) TR T ANGULAR BROAD CRESTED WEIR,”/ 

5 I0X, "<4)0UTCHER MOVABLE WEIR”//) 

GO TO 94 

4L I T ( Q1 MAX . LT . £ . 1 ) GO TO SO 

WRITE £6, 77) 

IF< IANS.EQ. 1 ) GO TO 303 
GO TO 94 

303 WRITE ( 3£, 77 ) 

7? FORMAT ( 1 OX, "THE RELEVANT STRUCTURE IS"/ 

I f OX, -RACIAL OR TAINTER GATE"//) 

GO TO 94 

50 I F ( ISED . EG . 3 ) GO TO 52 

IFIQIHAX.LT. 0.4JSO TO S4 
WRITE I 6, 7$ > 

WRITE < 32 , 79 ) 

I F f JANS. EG. 1 ) GO TO 304 
GO TO 94 

304 WRITE £32, 79) 

79 FORMATU OX , "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, ** £ 1 JRAD1 AL OR TAINTER GATE"/ 

5 JOX, ”<£)CRUMP OE GRUYTER ADJUSTABLE ORFICE”/ 

1 1 OX, *' £3>METER GATE"//) 

GO TO 94 

52 WRITER 6, e n 

IF£ IANS .EQ . 1 ) GO TO 30S 
GO TO 94 

305 WRITE £ 32, S 1 ) 

8 1 FORMATdOX, "THE RELEVANT STRUCTURES ARE"/ 

1 1QX,*‘£1 JRADIAL OR TAINTER GATE* 1 / 

1 1 OX, ” ( £ )RECT ANGULAR SHARP EDGED ORIFICE"//) 

GO TO 94 

S4 WRITE £ 6 , S3 ) 

1 F ( IANS.EQ. 1 ) GO TO 306 
GO TO 94 

306 WRITE! 3£ , 81 ) 

03 FORMAT £ 1 OX, "THE RELEVANT STRUCTURES ARE**/ 

1 1 OX, "(1 JRAOIAL OR TAINTER GATE”/ 

1 1 OX, “£2>CRUMP DE GRUYTER ADJUSTABLE ORIFICE**/ 

1 10X, “ ( 3 )METER GATE”/ 

1 10X, *‘£4)NEYRPIC MODULE “/ 

5 1 OX , " ( 5 ) RECTANGULAR SHARP EDGED ORIFICE"//) 

GO TO 94 

94 RETURN 

END 


SUBROUTINE TRPS £ Hi MAX , H2MAX, 01 MAX , 01 Ml N , P, B, ISED, IANS) 
£ F ( HI MAX . LT . £ . 0 ) GO TO 156 
IF HANS. EG. 1 ) GO TO 307 
GO TO 390 

07 WRITE (3£,B5) 

9 0 URITE£&,85> 

85 FORMATU OX, "THE RELEVANT STRUCTURE IS"/ 

i 1 OX, "CYLINDRI CAL CRESTED WEIR**//) 

GO TO 96 
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I*'. CONTINUE 

CO TO 56 

5^ IF * HI MAX . LT . t . 0 >G0 TO 58 

IFHANS.EO. J )C0 TO 308 
GO TO 391 

308 WRITE <32 ,87) 

391 UR I TE ( 6 , 87 ) 

87 FORMAT* TOX, "THE RELEVANT STRUCTURES ARE’’/ 

1 1 OX, “ <1 )CYLJNDRICAL CRESTED UEIRV 

3 1 OX , " * 2 ) TR l ANGULAR BROAD CRESTED WEIR"/ 

1 1 OX, “(3). ROUND NOSE BROAD CRESTED WET RV/) 

GO TO 96 

58 AML*H£'MAX/H J MAX 

WRITE**,*) ’AML*' , AML 

JF { ( AML . GT . 0 . 3 ) .AND. {AML.LT. 0 . S ) ) GO TO 60 
IF < (AML .GT . 0 .5) . AND . { AML.LT. 0 . 6 ) >GO TO 62 
IF (CAML.&T.O.S) .AND. (AML.LT. 0.7) > GO TO 64 
IF ( ( AML.GT . 0 . 7 ) . AND . (AML .LT . 0 .8) )GO TO 66 
IF << AML. GT. 0. S ) .AND. ( AML.LT. 0.95) ) GO TO 68 

CO 3 F ( I ANS . EG . 1 ) GO TO 309 

GO TO 392 

30 8 WRITE <38,89 ) 

3 VC UR I TEC 6, 89) 

88 FORMAT ( 1 OX, “THE RELEVANT STRUCTURES ARE*'/ 

3 1 OX , "ROUND NOSE HORIZONTAL BROAD CRESTED WEIR"/ 

1 1 OX, "THE RAMI JIN MOVABLE MEASURING OR REGULATING WE TR "/ 

1 J OX, ”TR (ANGULAR BROAD CRESTED WEIR"/ 

1 10X, “BUTCHER MOVABLE UE IR “ / 1 OX , " D I V I SOR " // ) 

GO TO 96 

6L I F * IANS. EG. 1 ) GO TO 35 0 

GO 70 393 

310 WRITE* 3£, 91} 

393 WRITE (6,91 I 

91 FORMAT (1 OX, "THE RELEVANT STRUCTURES ARE"/ 

1 1 0X, "ROUND NOSE HORIZONTAL BROAD CRESTED WEIR"/ 

3 T OX , "TR I ANGULAR BROAD CRESTED WEIR"/ 

1 1 OX, “BUTCHER MOVABLE WE IR " /I OX , “ D I V ISOR " // ) 

GO TO 96 

64 CFUANS.EG. 1 )GO TO 331 

GO TO 394 

311 WR I TE * 32 , 93 ) 

394 WR ITE( 6 , 93 ) 

93 FORMAT*? OX, "THE RELEVANT STRUCTURES ARE”/ 

1 1 OX, "ROUND NOSE HORIZONTAL BROAD CRESTED WEIR"/ 

3 1 OX, "TRIANGULAR BROAD CRESTED WEIR”/ 

1 1 OX, "BUTCHER MOVABLE WEIR"//) 

GO TO 96 

66 IF U ANS . EG . 1 )GO TO 312 

Gii TO 395 

31 £ WRITE <32, 95) 

395 UR ITE ( 6 , 95 ) 

9 5 FORMAT* 1 OX, "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, “ROUND NOSE HORIZONTAL BROAD CRESTED WEIR“/ 

3 3 OX , "TRIANGULAR BROAD CRESTED WEIR"//) 

GO TO 96 
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t.r 3 r ( ians-Cq. i >go to 313 

GO TO 396 

313 WRITE(3£, 97) 

396 WRITE (6, 97) 

97 FORMAT! 1 OX, "THE RELEVANT STRUCTURE IS“/ 

? 3 OX , "ROUNO NOSE HORIZONTAL BROAD CRESTEO WEIR"//) 

GO TO 96 
96 RETURN 

END 


SUBROUTINE SDNO < HI MAX, H2MAX, HI MI N, Q1 MAX, G1 MI N, GAMA, IANS) 
WRITE<*,*)'GAMA=' , GAMA 
3 F ( G1 MAX . GE .0.4) GO TO 1111 
3 FfQIMAX.LT. 0.4) GO TO ££££ 

1111 CALL SDRS i HI MAX , H2MAX, H1HIN.Q1 MAX, Q1MIN,P,B, IANS) 

GO TO 98 

aaaa i f <oama.lt. 56) go to 9 

WRITE! 6,33 ) 

IF( IANS.EQ. 1 ) GO TO 314 
GO TO 9 a 

314 WRITE (32, 31) 

31 FORMAT! I OX, "THE RELEVANT STRUCTURES ARE”/ 

1 1 OX, “TR3 ANGULAR SHARP CRESTED WEIR"/) 

GO TO 98 

9 UN=< ALOG! Q1 MAX ) -ALOG ( Q 1 M IN > >/! ALOG < HI MAX ) -ALOG ! H ! M IN > ) 

WRITE!*,*) ' UN= ' , UN 


IF 

£ (UN 

. GT. 0 . £5 ) 

.AND. 

(UN 

. LE . 

0 .7S) ) 

GO 

TO 

1 1 

) F 

£ (UN 

.GT. 0.7S) 

. AND. 

£ UN 

.LE. 

1 . £S ) ) 

GO 

TO 

12 

IF 

£ £ UN 

.GT. 1 ,£S) 

. AND . 

(UN 

.LE . 

1 . 7S > ) 

GO 

TO 

13 

3F 

£ £ UN 

. GT . 1 . 7S > 

. AND . 

£ UN 

.LE . 

2. £5) ) 

GO 

TO 

1 4 

IF 

£ UN . 

GT. £ . £S 3 

GO TO 

IS 







11 WR I TE C £ , 33 ) 

IFUANS.EQ. 1 ) GO TO 31S 
GO TO 98 

15 WR I TE ( 32 , 33 ) 

3 FORMAT ( 1 OX , "THE RELEVANT STRUCTURES ARE"/ 

1 )0X, ,! ( 1 JCIRCULAR SHARP EDGED ORIFIC”/ 

1 1 OX, “ ( £ ) RECTANGULAR SHARP EDGED ORIFIC"/ 

1 1 OX , ” ( 3 ) DANA I DEAN TUB"//) 

GO TO 9S 

£ WE ITE ( 6 , 3S ) 

I F ! IANS.EQ. 1 ) GO TO 316 
GO TO 98 

16 WRITER 32, 35) 

5 FORMAT! 1 OX, “THE RELEVANT STRUCTURE IS"/ 

1 1 OX, "PROPORT IONAL WEIR”//) 

GO TO 90 

13 UR ITE (6, 37) 

IF (IANS.EQ. I ) SO TO 317 
GO TO 98 

317 WRITE (3£, 37) 

37 FORMAT ( 1 OX, “THE RELEVANT STRUCTURE IS"/ 

1 1 OX, "RECTANGULAR SHARP CRESTED WIER’V/) 

GO TO 90 

14 UR ITE ( 6 , 39 ) 
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f F ( IANS. EG. J )C0 TO 318 
GO TO 98 

3 1 8 WRITER 32, 39 ) 

39 FORMAT < 1 OX, "THE RELEVANT STRUCTURE ISV 

1 J OX, "CIRCULAR WIER"//) 

GO TO 90 

15 URITEC 6, 41 ) 

f F < IANS. EG. 1 )G0 TO 319 

go to se- 
ns UR I TE ( 32 , 4 1 ) 

41 FORMATC 1 OX, "THE RELEVANT STRUCTURE IS'V 

J 1 OX, "TRIANGULAR SHARP CRESTED UEIR’ 1 //} 

GO TO 90 
90 RETURN 

END 

C 

SUBROUTINE SOPS { HI HA X, HEM AX, HI HI N, 01 MAX . 01 MIN, P, B , IANS ) 

GPN= (GIHAXJ/e 
Wfi I TE < * , * ) ' GRN® ' , QPN 
TFCQPN.LT.6.0) GO TO 16 
UR I TEC 6, 43 ) 

41 FORMAT C 1 OX, "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, "LARGE CAPACITY SPILLWAYS"/ 

1 1 OX, "(DUES STANDARD UEIR”/ 

1 1 OX , " ( 2 ) CYL I NDR I CAL CRESTED UEIR"/) 

GO TO 86 

16 CALL SU8CCH1MAX, HEM AX ,Q1HAX,G1MIN,P,E, IANS) 

GO TO 321 

86 I F { I ANS . EQ . 1 ) GO TO 3£0 

GO TO 321 

320 UR I TE C 32 , 43 ) 

321 RETURN 
END 

C 

SUBROUTINE SUBCC H ? MAX , HEMAX , G I MAX , G 1 M I N , P , B , I ANS > 

GAMA-Q1 MAX/G1 MIN 
URITEC*, *) ' GAMA= ' , GAMA 
IFCGAMA.GT.3S0) GO TO 17 
IFIP.LT. 0. 01 ) GO TO 19 

CALL SUSP 1 < H ! MAX , HcMAX , G 1 MAX , Q 1 M IN , P , B , I ANS ) 

GO TO S3 

17 I FCP .LT . G . 01 ) GO TO 18 

CALL SUBP2CH1MAX, H2MAX , G1MAX,G1M!N,P,B, IANS) 

GO TO 88 

IS IFC CQ1MAX .GT . 3 . 3) . OR . (HI MAX . GT . 1 .8) )G0 TO £0 

UR I TE ( & , 4S ) 

ire IANS .EG .1 ) GQ TO 322 
GO TO GO 

322 UR I TEC 32, 45) 

4 C . FORMAT U OX, "THE RELEVANT STRUCTURE IS"/ 

1 1 OX, “H- FLUME"//) 

GO TO 88 : 

59 IF < G1MAX.LT. 3. 3) GO TO 21 

UR I TE ( 6 , 47 ) 

I F C I ANS . EG . 1 ) GO TO 323 
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47 
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324 
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326 

4 9 


SO 

C 


327 

SI 


£.‘4 


fcO TO 6ft 
UR I TC l 32 , 47 ) 

FORMAT < J 0X , "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, - ( DPARSHAL FLUME ** /I OX , “ ( 2 ) LONG THROATED FLUME**//) 

CO TO CO- 

CALL SU8F2 ( H 1 MAX , HE MAX , Q ! MAX , Q J M IN , P , 6 , IANS ) 

GO TO 80 

I FCH1 MAX . LT . 1 . 36 ) GO TO 22 
UR I T£ ( £ , 147) 

IFUANS. EG. 1 ) CO TO 324 
CO TO 88 
URITE<3£, 147) 

FORMAT C 1 OX, “THE RELEVANT STRUCTURES ARE"/ 

1 1 OX , " * 1 ) MARSHAL FLUME "/ T OX , ”<£) LONG THROATED FLUME"//) 

GO TO 38 
AML=H2MAX/H1 MAX 
UR ITE C * , * ) 1 AML= ' , AML 
IF* AML . LT . 0. £S) GO TO £3 
WRITE* 6 ,£47) 

JFCIANS.EG. 1 ) GO TO 325 

GO TO S3 

UR I TEC 32.. 24?) 

FORMAT * 1 OX , "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, *' ( 1 MARSHAL FLUME "/I OX, *’ C2)L0NG THROATED FLUME**//) 

CO TO 68 
URI TE < 6 , 49 ) 

1FC IANS .EG! . 1 ) GO TO 326 
GO TO as 
UR I TEC 32, 49 ) 

FORMAT* 1 OX, "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, ** C 1 )WE)R SILL WITH RECTANGULAR CONTROLSECTION'*/ 

1 1 OX , " ( 3 JPARSHAL FLUME ”/ 1 OX , “ ( 4 ) H-FLUME ■// ) 

GO TO 3S 

RETURN 

END 


SUBROUTINE SUBR1 C HI MAX , H2MAX , Qt MAX, 01 MIN, P , B, IANS) 
GAMA*G 1 MAX/Q 1 M i N 

( F < * H2MAX/H1 MAX ) . CT . 0 . 66 ) GO TO £4 
IF*H1MAX. GT.£. OJGO TO £5 
IF*Q1MIN.LT. 0. 006 ) GO TO £6 
I F * GAMA . GT . 3£ ) £0 TO £7 
WRITE * 6 , 51 ) 

3FC IANS .EG. 1 >G0 TO 327 
GO TO o c 
WRITE C 32 , SI ) 

FORMAT* I OX, "THE RELEVANT STRUCTURES ARE"/ 

1 1 OX, “ C 1 ) RECTANGULAR CONTROL SECTION WITH"/ 

5 i OX , " * £ ) ROUND NOSE HORIZONTAL BROAD CRESTED WEIR"/ 
1 1 C-X, " C3)BR0AD CRESTED RECTANGULAR PROFILE UEIRV 

1 1 OX, ” * 4)FAIYUM WEIR"/ 

1 1 OX , ** C S ) TR 1 ANGULAR PROFILE TUO DIMENSION UEIR**/ 

1 1 OX, " (S)CYLINDRICAL CRESTED WEIR"/ 

1 1 OX, “ C7 )LONG THROATED FLUME**//) 

GO TO 8£ 

IF* H1MAX.LT. £.6) GO TO ! £5 
WRITE* 6, 52) 
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370 

3L 

1 

1 

1 •'.•r 

* C 


328 
1 S3 

1 

1 

1 

1 

126 


329 
1 33 

1 

? 


330 

ss 

1 

1 

1 

1 

26 


331 

137 

i 

1 

1 

t 


1 


I F ( JANS CO . 1 )CO TO 370 
GO TO 62 
WRITE ( 32 , S2 ) 

FORMAT ( 1 0X , "THE RELEVANT STRUCTURES ARE"/ 

1 OX , " ! 1 ) TR I ANGULAR PROFILE TWO DIMENSION UEIR"/ 

1 OX, "(2)L0NG THROATED FLUME"//) 

00 TO 82 

]F!G1M1N.LT.O.O06>GO TO 126 
URITE!6, 1S3) 

IF< IANS. EG. J )C0 TO 328 
GO TO 62 
WRITE (32, 133) 

FORMAT! 1 OX, “THE RELEVANT STRUCTURES ARE"/ 
10X,"(1)ROUN0 NOSE HORIZONTAL BROAD CRESTED WEIR”/ 

1 OX, ** !2)TRI ANGULAR BROAD CRESTED WEIR “/ 

1 OX, "(3)TR (ANGULAR PROFILE TWO DIMENSION WEIR"/ 

1 OX , ** < 4 ) LONG THROATED FLUME**//) 

GO TO 82 
WR I TE < 6 , 133) 

I F < IANS .EG . 1 > GO TO 329 

GO TO 8 2 

UR I TE { 32, 1 S3 > 

FORMAT ( 1 OX , "THE RELEVANT STRUCTURES ARE”/ 

1 OX, “ !£)CYLINDRICAL CRESTED UEIR"/ 

1 OX , * ( 4 ) LONG THROATED FLUME"//) 

GO TO 82 
UR I TE ( 6 , SS J 

I F ( IANS. EG. ? ) GO TO 330 
GO TO 82 
UR I TE ( 32 , 33 > 

FORMAT ! 1 OX, “THE RELEVANT STRUCTURES ARE**/ 

1 OX, * (1 )TR I ANGULAR PROFILE FLAT-V-UEIR”/ 

1 OX , " ! 2 ) LONG THROATED FLUME"/ 

1 OX, ” ( 3)CYL INDRICAL CRESTED UEIR"/ 

1 OX , " ( 4 > LONG THROATED FLUME"/) 

GO TO 82 
WRITE! 6 , 137) 

IF!) ANS .EQ . 1 )G0 TO 331 
GO TO 32 
WRITE! 32 , 1 S7 > 

FORHATUOX, "THE RELEVANT STRUCTURES ARE"/ 

1 OX, " t 1 ) TRIANGULAR BROAD CRESTED UEIR"/ 

IOX,"(£)FAI YUM WEIR”/ 

1 OX , “ I 3 > V-NOTCH UEIR SILL"/ 

1 OX , " < 4 } LONG THROATED FLUME’’//) 

GO TO 82 
WRITE! 6, S7) 

! F < IANS.EO. 1 >G0 TO 332 
GO TO 82 
UR I TE < 32 , 37 ) 

FORMAT! 1 OX, “THE RELEVANT STRUCTURES ARE"/ 

5 OX ,’’( 1 )TR (ANGULAR BROAD CRESTED WEIR’’/ 

1 OX , “ ! £ > CYLINDR I CAL CRESTED UEIR"/ 


i 



ra uj r> r:« uj 
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-4 4*. 1*1 u 


1 IPX, ” < 3 ) F A I YUM WEIR"/ 

1 10X, " ! 4 ) TR 1 ANGULAR PROFILE FLAT-V-UEIR"/ 

J ! OX , " ! 5 ) LONG THROATED FLUME”/) 

GO TO 32 
RETURN 
END 

OUBKOUT I NE SUB PE ( H J MAX , HEM AX, G ? MAX , QIMJN, F , B , IAN S~ 
I F( H2MAX/H1 MAX . LT . 0 . 66 >GG TO 22-4 
WRITE! 6 , 2SS ) 

IF! IAN5.EQ. 1 )GO TO 333 
CO TO £92 
WRITE(32, £5?) 

FORMAT! 1 OX, "THE RELEVANT STRUCTURE IS"/ 

5 1 OX, ”TR IANCULAR BROAD CRESTED WEIR"//) 

GO TO 292 

IF! H1MAX.LT. 2. 0>GO TO 226 
UR I TE i 6 , 257 5 
!F< IANS. EG. 1 >GO TO 334 
GO TO £92 
WRITE! 32, 2S7) 

FORMAT! 1 OX, “THE RELEVANT STRUCTURES ARE*'/ 

? 1 OX , ”TR ! ANGULAR PROFILE FLAT-V-WE IR”/ 

1 1 OX, “CYLINDRICAL CRESTED WEIR*'//) 

GO TO 292 
£26 IFIQIMFN.GT. 0. 006)G0 TO 228 

WR I TE ( 6 , 259 ) 

IF! IANS .EG . 1 ) GO TO 33S 
GO TO 292 

335 UR I TE ( 32, 259 > 

£59 FORMAT!) OX, "TRIANGULAR BROAD CRESTED WEIR"/ 

1 1 OX, “V-NOTCH UEIR SILL"//) 

GO TO 292 

£28 WRITE <6, 26)) 

IF! IANS .EG!. 1 ) GO TO 336 
GO TO 292 

336 UR I TE ! 32 , 26 1 ) 

26) FORMAT! 1 OX, "THE RELEVANT STRUCTURES ARE”/ 

1 1 OX , “ £ 1 ) TR I ANGULAR BROAD CRESTED UEIR"/ 

) ) OX," ! 2) V-NOTCH UEIR SILL”/ 

1 i 0 X , “ £ 3 > C YL I NDR I CAL CRESTED UEIR"//) 

CO TO 29c 
292 RETURN 

END 
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SELECTED RELEVANT MEASURING DEVICES BY 
PROGRAM SELECT FOR A GIVEN INPUT. 

EXAMPLE 1 

G 1 MAX THE MAXIMUM DISCHARGE IN CUBIC METER S/SEC 
GTMIN THE MINIMUM DISCHARGE IN CUBIC METERS/SEC 
Y 1 MAX THE MAXIMUM DEPTH OF WATER IN METERS 
Y1MIN THE MINIMUM DEPTH OF WATER IN METERS 
P HEIGHT OF CREST IN METERS 
DHMX MAXIMUM HEAD LOSS IN METERS 

B BREATH OF BOTTOM OF CONTROL SECTION IN METERS 
INPUT DATA 


( 1 ) Y1 MAX 

SB 

3 . 00000 

METERS 

CE ) 

Y1 MIN 

= 

£ . 40000 

METERS 

(3) 

Q 1 MAX 

* 

3 . 00000 

CUBIC 

( 4 ) 

Q1 MIN 

St 

. 03000 

CUBIC 

(5) 

P 

3= 

1 . 40000 

METERS 

( 6 ) 

B 

ss 

E . 00000 

METERS 

C 7 ) 

DHMX 

= 

. 60000 

METERS 


METERS/SEC 

METERS/SEC 


DISCHARGE MEASURING STRUCTRES 
SEDIMENT OR TRASH TO BE PASSED 
TRASH NOT TO BE PASSED 
SUBMERGED FLOW CONDITION 


THE RELEVANT STRUCTURES ARE 
C 1 JTRIANGULAR BROAD CRESTED WEIR 
CSJCYLINDRI CAL CRESTED WEIR 
C 3 ) FA I YUM WEIR 

( 4 > TR I ANGULAR PROFILE FLAT-V-WEIR 
t S ) LONG THROATED FLUME 


EXAMPLE £ 

Q1 MAX THE MAXIMUM DISCHARGE IN, CUBIC METERS/SEC 

Q1MIN THE MINIMUM DISCHARGE IN CUBIC METERS/SEC 

Y 1 MAX THE MAXIMUM DEPTH OF WATER IN METERS 

Y1MIN THE MINIMUM DEPTH OF WATER IN METERS 

P HEIGHT OF CREST IN METERS 

DHMX MAXIMUM HEAD LOSS IN METERS 

B BREATH OF BOTTOM OF CONTROL SECTION IN METERS 


INPUT DATA 


( 1 ) Yt MAX = 
(£) YtMIN = 

(3) Q 1 MAX = 

(4) Q1MIN = 

< 5 ) P = 

( 6 ) B * 

(7) DHMX * 


2.00000 METERS 
1 .50000 METERS 
.30000 CUBIC METERS/SEC 
.03000 CUBIC METERS/SEC 
1 .00000 METERS 
1 .60000 METERS 
.50000 METERS 



discharge measuring structres 

SEDIMENT OR TRASH NOT TO BE PASSED 
TRASH NOT TO BE PASSED 
SUBMERGED PLOW CONDITION 


THE RELEVANT STRUCTURE IS 
TRIANGULAR SHARP CRESTED WEIR 


EXAMPLE 3 

Q 1 MAX THE MAXIMUM DISCHARGE IN CUBIC METERS/SEC 
Q1MIN THE MINIMUM DISCHARGE IN CUBIC METERS/SEC 
Y 1 MAX THE MAXIMUM DEPTH OF WATER IN METERS 
Y1MIN THE MINIMUM DEPTH OF WATER IN METERS 
P HEIGHT OF CREST IN METERS 
DHMX MAXIMUM HEAD LOSS IN METERS 

B BREATH OF BOTTOM OF CONTROL SECTION IN METERS 


INPUT DATA 


< 1 >Y1 MAX = 

(C) r j ni n = 
(3) 0 1 MAX « 

<A) CUMIN = 

< 5 > P * 

C 6 J B ■= 

( 7 > DHMX ■ 
<B3PLX= .so 


£.00000 METERS 
1 .40000 METERS 

6.00000 CUBIC METERS/SEC 
.30000 CUBIC METERS/SEC 

1 . E0000 METERS 

2 . 00000 METERS 
.20000 METERS 


DISCHARGE REGULATING STRUCTURES 
SEDIMENT OR TRASH NOT TO BE PASSED 
TRASH NOT TO BE PASSED 
SUBMERGED PLOW CONDITION 


THE RELEVANT STRUCTURES ARE 

(I) MOVABLE ROUND NOSE HORIZONTAL BROAD CRESTED WEIR, 
CS) MOVABLE TRIANGULAR BROAD CRESTED WEIR, 

(3) BUTCHER MOVABLE WEIR 
tA> DIVISORS 


EXAMPLE A 

Q 1 MAX THE MAXIMUM DISCHARGE IN CUBIC METERS/SEC 
CUMIN THE MINIMUM DISCHARGE IN CUBIC METERS/SEC 
Y 1 MAX THE MAXIMUM DEPTH OF UATER IN METERS 
Y1M1N THE MINIMUM DEPTH OF WATER IN METERS 
P HEIGHT OF CREST IN METERS 
DHMX MAXIMUM HEAD LOSS IN METERS 

B BREATH OF BOTTOM OF CONTROL SECTION IN METERS 


i 15 



INPUT DATA 


( 1 ) Y I MAX « 

( Z ) Y 1MIN « 
(3) (UMAX 
< A ) a IMIN * 
(3> r * 

<61 B « 

< T ) DMRX * 

< 8 JTLX* I . ZO 


1 .£0000 METERS 
.80000 METERS 
C. 00000 CUBIC HETER/3EC 
.40000 CUBIC METEM/'SEC 
.30000 METERS 
Z . 00000 METERS 
.30000 METERS 


DISCHARGE REGULATING STRUCTURES 
SEDIMENT OR TRASH NOT TO BE TASSED 
TRASH NOT TO BE MASSED 
SUBMERGED TLOW CONDITION 


THE RELEVANT STRUCTURES ARE 
< 1 JRADIAL OR TAINTER GATE 
ICJCRUnr DE GRUTTER ADJUSTABLE ORriCE 
< 3 I METER GATE 


EXAMPLE 5 

0) 1 MAX THE MAXIMUM DISCHARGE IN CUBIC METERS/SEC 
CUMIN THE MINIMUM DISCHARGE IN CUBIC METER5/3EC 
Y 1 MAX THE MAXIMUM DEPTH OF WATER IN METERS 
Y1MIN THE MINIMUM DErTH OF WATER IN METERS 
r HEIGHT OF CREST IN METERS 
DHMX MAXIMUM HEAD. LOSS IN METERS 

B BREATH OF BOTTOM OF CONTROL SECTION IN METERS 
INPUT DATA 


< t ) Y 1 MAX = 

( £ ) YT MIN » 

< 3 ) 01 MAX * 

( 4 ) Q! MIN = 

(5) P * 

(6) B = 

(?) DHMX = 

< 8 ) FLX= 1 .£0 


1 . 1 0000 METERS 
.80000 METERS 
.30000 CUBIC METERS/SEC 
.03000 CUBIC METERS/SEC 
.70000 METERS 
£ . 00000 METERS 
.40000 METERS 


DISCHARGE REGULATING STRUCTURES 
SEDIMENT OR TRASH NOT TO BE PASSED 
TRASH NOT TO BE PASSED 
FREE FLOW CONDITION 


THE RELEVANT STRUCTURES ARE 
(l)RADIAL OR TAINTER GATE 
< £ ) CRUMP DE GRUYTER ADJUSTABLE ORIFICE 

( 3 ) METER GATE 

( 4 ) NEYRP I C MODULE 

( 5 ) RECT ANGULAR SHARP EDGED ORIFICE 
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C PRGRAM DESIGN 

c a rRoenfln rcm design or weirs 

REAL K 
INTEGER CN 

OrCNtgNIT«CD / riLE*'DC3IGN.OUT’ ) 

WRITE!©, t 13) 

113 FORMAT! 1 OX, 'DESIGN OF WEIRS’) 

WRITE! 6, 115) 

115 FORMAT! 1 OX, 'LIST OF WEIRS WITH CODE NUMBER'/ 

1 JdX,'l LINEAR PROPORTIONAL WEIR'/ 

1 10)|, '2 • QUADRATIC WEIR'/ 

1 1 OX , ' 3 LOGRITHMIC WEIR'/ 

1 1 OX, '4 EXPONENTIAL WEIR'/ 

I 1 OX , '5 NEW EA5ELES5 WEIR NBWE ' / 

1 1 OX , '6 WES STANDARD WEIR'/ 

7 1 OX , '7 ROUND NOSE BROAD CRESTED WEIR'/ 

7 1 OX , '8 ROM J I N MOVABLE WEIR'/ 

7 1 0X , ' 9 TRIANGULAR BROAD CRESTED WEIR'/ 

1 10X, '10 BROAD CRESTED RECTANGULAR WEIR'/) 

WRITE(6, 7 1 7 J 

117 FORMAT! 7 OX, 'ENTER THE CODE NUMBER DESIRED') 

READ ( 6 , * ) CN 
WRITEC6, EEC) 

!S£ FORMAT! 7 OX, ' G 1 MAX=MAXI HUM DISCHARGE IN CUBIC METERS/SEC'/ 

7 1 OX, 'Q1MIN=MINIMUM DISCHARGE IN CUBIC METERS/SEC'/ 

7 1 OX , ' B=WIDTH OF WEIR IN METERS'/ 

1 1 OX, 'H1MAX-MAXIMUM HEAD IN METERS'/ 

1 1 OX, '^PROPORTIONALITY PARAMETER'/ 

1 1 OX, 'FL=LENTGH ALONG FLOW DIRECTIONS'/ 

7 1 OX j ' YD=DOWN STREAM Y COORDINATE PROFILE'/ 

1 1 OX , ' OF WES-SPILLWAY'/ 

1 70X,'YU=D0WN STREAM Y COORDINATE PROFILE'/ 

1 1 OX , ' OF WES-SPILLWAY'/ 

1 7 OX , ' X=X COORDINATE OF PROPORTIONAL WEIR'/ 

1 1 OX , ' Y“ Y COORDINATE OF PROPORTIONAL WEIR'/) 

WRITE (*,310) 

10 FORMAT! 1 OX, "DO YOU WANT TO STORE RESULTS IN A FILE NAMED", 

1 " DESIGN. OUTVIOX, "IF YES WRITE t OTHERWISE 2") 

READ! *, *>IANS 

IF! IANS.EQ. 1 ) GO TO 375 

GO TO 376 

75 UR I TE ( £0 , £22 ) 

76 IF ! CN . EQ . 1 > GO TO 119 
IF ! CN . EQ . £ ) GO TO 1£1 
IF! CN . EQ . 3 ) GO TO 1£3 
IF ! CN . EQ . 4 ) GO TO 1 £5 
I F ( CN . EQ . 5 ) GO TO 1£7 
IF ! CN . EQ . 6 ) GO TO 1£9 
IF ( CN . EQ . 7 ) GO TO 131 
IF ! CN . EQ . 8 ) GO TO 133 
IF ! CN . EQ . 9 ) GO TO 135 
I F ! CN . EQ .10) GO TO 137 
CALL SUTRO! IANS) 

GO TO 100 


19 



» 10 

:1 CALL GUADRUANS) 

CO TO 100 

'3 CALL LOCK IANS) 

CO TO 100 

5 CALL EXPO(IANS) 

CO TO 100 

'.7 CALL NEWS ( IANS) 

CO TO 100 

:s CALL WES < IANS ) 

CO TO 100 

(1 CALL RNBU(IANS) 

GO TO 100 

!3 CALL RMW(IANS) 

GO TO 100 

t5 CALL TBU(IANS) 

GO TO 100 

57 CALL BCR! IANS) 

GO TO 1 00 

10 STOP 

END 


SUBROUTINE SUTRO(IANS) 

REAL K 

WRITE 56,100) 

10 FORMAT ! 1 OX, 'WELCOME TO THE DESIGN OF'/ 

1 10X, 'LINEAR PROPORTIONAL WEIR'/) 

PAUSE 

WRITE { 0 , 11 1 > 

1 FORMAT! 1 OX, 'ENTER THE VALUE Q 1 MAX , Q 1 M I N, Y1 MAX , B , DY ' ) 

READ ! 6 , * ) Q 1 MAX , Q 1 MI N , Y 1 MAX , B , DY 

G*9.81 

CD=0.63 

A=!GHMIN/!CD*2./3.*!E.*G)**0.5*B) )**CE./3. ) 

WRITE! *, * ) A 
K«OD*B*!£. 0*G*A ) **Q . 5 
WRITE! * , * )K 
H1MAX*Q1 MAX/K+A/3 . 

P*T1 MAX-H1 MAX 

WRITEI6,*) INPUT DATA ' 

WRITE! 6 , 1 33)QIMAX,01MIN,B, Y1MAX, CD, G 
I 33 FORMAT! IX, ' Q 1 MAX* ',F10.4,EX,'Q1MIN=',F10.4,EX,'B=',F10.4, 
1 2x , ' Y1MAX=' , FI 0 . 4 , EX , 'CD=' , F 1 0 . 4 , EX , ' G= ' , F 1 0 . 4/ ) 

IF! IANS .EQ . 1 ) GO TO 31 1 
GO TO 31 E 

11 WRITEIEO, * ) ' INPUT DATA ' 

WRITE! EO, 1 33 ) Q 1 MAX , Q 1 MI N , B , Y 1 MAX , CD , G 

UR I TE ! £0 , * ) ' OUTPUT FOLLOWS ' 

WRITEIEO, £00)A,K,H1MAX,P 

IE WRITE!6,*) ' OUTPUT FOLLOWS ' 

WRITE!6,£00)A,K,H1 MAX , P 

0 0 FORMAT! IX, ' A® ' , FI 0 . 4 , EX , 'K* ' , FI 2 . S, EX, ' HI MAX* ' , F 1 2 . 5 

1 , EX , ' P— ' F 1 0 .4/) 

Y»=0.0 

0 Y*=Y + DY 

X = 0 . 5*B* ! 1 . 0 -E . 0 / 3.141 59 £7* AT AN I SORT I Y/A ) ) ) 

IF ( Y . GT . HI MAX ) GO TO EO 
UR I TE ! 6 , 1 05 ) Y,X 



I i 9 

Ifi IANS 10 1 >00 tO 31 3 
CO TO 31 « 

13 write tec , i os ) t , x 

os f ormat 1 5 * . ’ t * ■ , e i e 5 , f x , ' x* * , r i r 5/ ) 

1 A CO TO 10 

o RtTUWN 

END 

SUBROUTINE QUAORt I ANS » 

WRITE t©,100) 

00 FORMAT ( 1 OX , * WELCOME rOR DE5ICN OF ' / 

1 1 ox, 'QUADRATIC WEIR '/) 

FAU5E 

WR I TE (6,111) 

11 FORMAT C 1 OX , ' ENTER THE VALUE Q 1 MAX , Q 1 M IN , T 1 MAX , B , DY ' ) 

READ ( 6 , * ) G I MAX , 0 1 W I N , T I MAX , B , DY 
G*9 .81 
CD=0 . 63 

A-CaiMIN/(E./3.*CD*(2.*G)**0.5*Bn**C2./3.> 

K^CO^B" (C . 0*G*A ? **-0 . 5 

FK= 1 . 1E4T*CD*B*A* ( ( E . 0*G ) **0 .5) 

H1MAX*CQ1MAX/FK)**CTC , *A/3, 0 
r*Tl MAX-HT MAX 
H 0*H 1 MAX-A 

WRITEC6, * ) ' INPUT DATA ' 

UR I TE < 6 , 133)Q1MAX,Q1MIN,B, Y1MAX, CD,C 

1 33 FORMAT ( IX, ' Q1 MAX- ' ,F1 0 , A, EX, ' Q 1 M I N- ' , F 1 0 . 4 , EX , 'B*' ,F1 0 .4 

1 ,Ex , ' Y 1 MAX** ' , FI 0 . 4 , BX, ' CD* ' , FI 0 . 4 ,EX, 'C*' ,F1 0 . 4/) 

I F ( IANS.EQ. 1 )CO TO 315 
GO TO 316 

15 WR I TE (20, * ) ' INPUT DATA ' 

WR I TE < £0 , 1 33 ) Q I MAX , 0 1 M I N , B , Y 1 MAX , CD , G 

WRITE ISO,*) ' OUTPUT FOLLOWS ' 

WRITE(£0,£00)A, FK, HI MAX, P 

16 UR I TE 16,*)' OUTPUT FOLLOWS ' 

URITE(6,E00)A,FK,H1MAX,P 

DO FORMAT ( 1X,'A*',F10.4,EX,'FK*',F1E.5,£X, 'HI MAX*' , FIE .5 

1 , 2X , ' P= ' F 1 0 . 4/ ) 

Y*Q . 0 

) Y = Y + D Y 

D*(1 0~£ 0/3 1 4 1 59S7*ATAN< SORT < Y/A) ) ) 

E = 6 . /3 . 14 1 S9E7*SGRT ( Y/A ) / ( 1 . +3 . * Y/A ) 

X* 0 5*B* ( D-E ) 

IF ( Y GT HI MAX ) GO TO EO 
WR I TE ( 6 , 1 OS ) Y . X 
I F < IANS EG 1 ) GO TO 317 
GO TO 3 1 e 

7 WRITE < EC , 1 05 ) Y,X 

FORMAT t EX , ' Y * ' , F 1 £ 5 , EX , ' X* ' , F 1 £ . 5/ ) 

GO TO 10 
RE TURN 
END 


3 


100 


1 1 1 


1 0 £ 


1 33 


31 3 


3£0 

£00 

1 

C 

30 


321 
1 0S 

322 
20 


SUBROUTINE LOCTC IANS* 

REAL K 

WRITE 16*100* 

FORMAT C 1 OX # f WELCOME TO THE DESIGN OF'/ 

1 OX* ' L0GR I THMI C WEIR f /1 
PAUSE 

WRITEC6* 1111 

FORMAT C 1 0X # * ENTER THE VALUE G 1 MAX > G 1 M I N , Y 1 MAX , B , DY ' 1 
READC6,*>G1 MAX , G 1 MIN, Y 1 MAX , B , DY 
G*«9. 81 
CD*0 . 63 

A*CQtMIN/CD*£ . /3. * (£ . *G ) ** . 5*B ) ** ( 2 . /3 . 1 

FFN* . 5 

FFN*0 .54-0 . 1 

CALL NRF<SG,FFN1 

WR 1 TE < * , * 1 SG 

K*G1 MIN/ALOGC SG ) 

WRITE( *, * )Q! MIN 

H0= (EXP (Q1 MAX/K )-SG)*FFN*A 

H1MAX*H0+A 

WR ITE( * , * )H1 MAX 

IFC HI MAX . LT . Y1 MAX ) GO TO 102 

r*rt MAX~*H1 MAX 

WRITCCS,*)' INrUT DATA * 

WR I TE ( 6 # 1 33)01 MAX , Q 1 M I N , B , Y 1 MAX , CD , G 

FORMAT C IX* r Q 1 MAX* ' , FI 0 . 4,2X, ' Q 1 M I N* ' , F 1 0 . 4 , 2X , ' B* ' , F 1 0 . 4 

* ex * ' Y1 MAX*' * FI 0 . 4* ex, ' CD* # , F 1 0 . 4 * 2X , ' G* ' , F 1 0 . 4/ ) 

IFC IANS. EG. 1 >G0 TO 319 

GO TO 320 

WR I TE (CO* * > * INrUT DATA 

WRITE (20, 1 33 1QTMAX, Q1MIN,B,Y1 MAX , CD * G 
WRITE (20* 200 1A,K,H1MAX,P 

WR I TE ( 6 , * 1 ' OUTPUT FOLLOWS ' 

WR I TE C6>£001A,K,H1 MAX * P 

FORMAT ( IX, ' A* ' , F 1 0 . 4 , 2X , ' K= ' * F 1 2 . 5 , 2X * 'HI MAX*' # FI 2. 5 

* 2X * 'P='F1 0.4/1 

Y = 0 0 
Y*T+DY 

E*1 .-C./3. 14*ATANCSQRT( Y/A1 1 
BF = £ . /3 . 1 4* ( SORT ( Y /A > > / ( SG*FFN+ Y/A 1 
BG=t C 5G^FFN+Y/A 1 ♦* . 5+ C Y/A1** .51 
BH=C CSG*FFN+Y/A1 **0 . 5-CY/A1** .51 
H*BG/BH 

AG=5G*FFN^C3 . 1 4^ CSG^FFNi'Y/'AJ 1**1 . 5*AL0GCH1 
X=0 .5*B*(E~BF-AG> 

IFCX . LT . 0 . 0 ) GO TO 20 

WR I TE C 6 > 1 05) Y * X 

IFC IANS . EG . 1 ) GO TO 321 

GO TO 322 

WR ITE ( 20 , 1 05 ) Y , X 

FORMAT (5X, 1 Y= f , F 1 2 . 5 , 2X * 1 X« ' * F 1 2 . 5/ ) 

GO TO 30 

RETURN 

END 



c 

370 

300 

C 

C 

1 00 

1 

1 1 1 


133 

1 

323 

3£4 

£00 

1 

40 

C 


SUBROUTINE NRF<3G, FN> 

SC-5.0 

FT- ( E . 0/3 . 0)/(SC*FN) 

FX«SG-EXP< FT > 

DFX-1 . 0+(FT/3G ) *EXP t FT ) 

5CI-SC -FX/DFX 

IF(ABS(5GI~5G ) . LE . 0. 0005)GO TO 300 

SG-SG 1 

CO TO 370 

RETURN 

END 


EXPONENTIAL WEIR 
SUBROUTINE EXFO(IANS) 

REAL K 

WRITE ( 6 , 100 ) 

FORMAT t 1 OX , 'WELCOME TO THE DESIGN OF'/ 

1 OX, 'EXPONENTIAL WEIR'/) 

PAUSE 

WRITE (6, 1 1 1 ) 

FORMAT! 1 OX, 'ENTER THE VALUE Q 1 MAX , Q 1 M I N , Y 1 MAX , B , CD , DY ' ) 
READ ( 6 , * )Q1 MAX, Q1 MIN, 71 MAX, B, CD, DT 
FLEM= 1 . 5 
G*9 .SI 

A* (OI MIN/CD*£. /3 . * (E . *G ) ** . 5*B )**(£. /3 . ) 
HO=A/FLEM*ALOG(G1MAX/GIMIN) 

HI MAX-HO+A 
P — Y 1 MAX-H 1 MAX 

WRITE C 6 , * ) ' INPUT DATA ' 

WRITE! 6, 133)01 MAX, Q1 MIN, B, Y1MAX,CD,G 

FORMAT ( IX, '01 MAX-' ,F10.4,£X, ' 0 I MI N* ' , F 1 0 . 4 , EX , ' B« ' , F 1 0 . 4 
,Ex, 'Y1 MAX-', FI 0.4, EX, ' CD- ' , FI 0 . 4, EX, 'G-', FI 0.4/) 


IF( IANS . EQ . 1 )G0 TO 3E3 
GO TO 324 

WR I TE ( £0 , * ) ' INPUT DATA 

WR I TE (20, 133 JQ1MAX, Q1MIN, B, Y1MAX, CD, G 

WRITEtEO,*)' OUTPUT FOLLOWS ' 

WRITE(£0,£00)A,H0,H1MAX,P 

UR I TE ( 6 , * ) ' OUTPUT FOLLOWS ' 


WRITE<6,£00)A,H0,H1MAX,P 

FORMAT < IX, 'A=' ,F1 0 .4,£X, ' HO- ' , FI £ . 5 , £X , 'HI MAX-' ,F1£.5 
,£X, ' P= ' F 1 0.4/) 

Y-0 . 0 
Y-Y+DY 

I F ( Y . GT . H 1 MAX ) GO TO 50 
OM-SGRT ( Y*FLEM ) 

CALL ERR (OM, TRAP) 


E=1 .-£./3.14*ATAN(SGRT(Y/A) ) 

H=£ . 0*01 MIN*FLEM** 1 . 5*EXP ( FLEM*Y ) / ( CD* ( £*9 . 81 )**.5>*TRAP 

X-0 . 5*B*E+H 

TYPE *, H, Y, OM, TRAP 

UR I TE ( 6 , 1 05 ) Y , X 

IF ( IANS . EQ . 1 )G0 TO 3£5 



CO TO 326 

325 WR1TE<20, 1 05) Y,X 

IDE FORMAT (5X, * Y-* ,F12 5,EX, * X* * , F 1 E . 5/ > 

326 60 TO 40 

50 RETURN 

END 

C 

SUBROUTINE ERR ( B , TRAP ) 

C PROGRAM TRPEZO I DAL 

C PROGRAM TO EVALUATE THE INTEGRAL OF EXP(X**-2.) 

C BETBEEN THE LIMITS A TO B 

C RULE OF INTEGRATION BASED ON N SUBINTERVALS 

C OR N+1 NODAL POINTS. THE VALUE OF A,B/ AND N 

C ARE TO BE READ AND THE ITEGRAND IS WRITTEN 

C AS A FUNCTION OFSUBPROGRAM 

C 

A* 0 . 0 
N* 1 7 

DXMB-AJ/N 
SUM = 0 . 0 
M*=N-1 

DO 20 1*1 , M 

X=A+I*DX 
SUM*SUM+r ( X ) 

20 CONTINUE 

TRAr*=0 . 5*DX*(F(A)+F(B)+2.*SUM) 

RETURN 

END 

C 

FUNCTION F< X ) 

F=1 . /EXP ( X**2 . ) 

RETURN 

END 

C 

C BASELESS WEIR C NBWE ) 

SUBROUTINE NBW2CIANS) 

REAL K 

WRITE (6,100) 

100 FORMAT ( 1 OX, 'WELCOME TO THE DESIGN OF'/ 

1 1 OX, 'BASELESS WEIR ( NBW£ ) ' / ) 

PAUSE 

WR I TE ( 6 ,111) 

111 FORMATC 1 OX, 'ENTER THE VALUE Q1 MAX , H 1 MAX , Y 1 MAX , T , CD , DY ' > 

READ ( 6 , * )G1 MAX , H 1 MAX , Y 1 MAX , T , CD , DY 
G=9 . 8 1 

AP=SQRTCH1MAX)*AL0G( 1 . +H1 MAX/T ) 

K=Q 1 MAX/ AP 

R = 3 . 1 41 5* CD* (SQRTtE . 0*G> )*T 
BU=K/R 

UR I TE ( 6 , * ) ' INPUT DATA ' 

URITE16, 133)Q1MAX,H1MAX,T,CD,DY 

FORMAT C IX, * Q 1 MAX- ' , FI 0.4, EX, ' H 1 MAX= ' , F 1 0 . 4 , EX , ' T* ' , FI 0 . 4 
1 ,£X, ' CD- ' , FI 0.4, EX, ' DY= ' , FI 0.4/) 
write( * , * )BW 
IF ( IANS . EQ . 1 ) GO TO 3E7 


1 33 



32 ?? 


3m 

c 

40 


C 


3E5 
1 05 
330 
50 

C 

C 

C 

C 

C 

C 

C 

C 

C 


£0 


C 


1 


c 

c 


100 

t 


CO TO 3£8 

WRIYEfEO,*? INPUT DATA — 

WR 2 T£ t £0 * 133)01 ItAX , H 1 ff AX , ? # CO * JDT 
WRZTEfCO,*) — OUTPUT FOLLOWS # 

UftlTE<6,*> # ' OUTPUT FOLLOWS--*- # 

Y*Y+DY 

XflY.GT.Hl MAX 1 GO TO 50 
CALL INTEGCY/T, TRAP > 

X*BU* TRAP “ 

WRITEC6,1D5> Y , X 
IFCIANS.EG. 1 >00 TO 3E9 

00 TO 330 
WRITECEO , 1 053 Y , X 

FORMAT CSX, 1 Y a * ,P1E .5# EX , r X* ' , F 1 £ . 5/ ) 

GO TO 40 

RETURN 

END 

SUBROUT INE I NTEG < Y~ T,TRAP > 

PROGRAM TRFEZOIDAL 

PROGRAM TO EVALUATE THE INTEGRAL 
BETBEEN THE LIMITS A TO B 

RULE OF INTEGRATION BASED ON N SUBINTERVALS 
OR N*M NODAL POINTS „ THE VALUE Or A,B, AND N 
ARE TO BE READ AND THE ITEGRAND IS WRITTEN 
AS A FUNCTION OFSUBFROGRAM 

A*-0~5 
B*0. 42999 
N*33 

DX a ( B-A ) /N 

5UM*0.0 

M*N~1 

DO £0 1*1 , M 
U= A+ I *DX 

SUM=SUM+FF (U,TiY) 

CONTINUE 

TRAP=0 .5*DX*(FF( A, T, Y)+FFCB, T, Y>+£. *SUM> 

RETURN 

END 

FUNCTION FF ( U7 T 7 77 

FF = SQRT i (U + 0 . 5 )/( 0 . 5-U) ) * ( 1 . 0/ ( < Y/T ) * ( U+0 .5)+1.0>**£.0+ 

0 . 5/ ( <Y/T)*(U+0.5)+t .0)) ( 

RETURN 

END 

WES- ST ANDARD SPILLWAY 
SUBROUTINE WES(IANS) 

WRITE (6,100) 

FORMAT ( 1 OX, * WELCOME TO THE DESIGN OF'/ 

1 OX, 'WES-STANDARD SPILLWAY '/> ; 



PAUSE 

WRJTECS* 101 ) 

i oi format* % ox, 'enter gimax # gimin # b, yimax # /i 

RE AD C 5 , * )Q1 WAX , Cl 1 HI H , B , Y 1 MAX 
WRITEC6, 38)G1MAX, GTMXN,B, Yf MAX 

38 FORMAT l ICX, * INPUT DATA. V 

1 1 OX, 1 MAXIMUM DISCHARGE G f MAX* ' , F 1 0 . 4 , £X , ' CUBIC METERS * / 
1 1 OX, -MINIMUM DISCHARGE 0 f H I N* ' , F 1 0 , 4 , EX , 'CUBIC METERS'/ 
1 10X, 'BREATH OF WEIR B* ' , F 1 0 . 4 , EX , ' METERS ' / 

1 1 OX * ' MAXIMUM DEPTH OF WATER Y T MAX* ' , F f 0 . 4 # EX , * METERS ' / ) 

WRITEC6, 1 O E) 

108 FORMAT ( 1 0X, 'ENTER TRIAL VALUE OF H 1 MAX ' / ) 

READ ( 5 , * ) H 1 M A X 
DO 111 1*1,10 

QBMAX*G1 MAX/B 
GBMIN*Q1 MIN/B 
VA*QBMAX/Y1 MAX 
HA*(VA**£. 0>/CE.0*9.81 ) 

EN*Y1 MAX+HA 

P*EN~HDMAX 

PHD«r/HDMAX 

C IF < < PHD . GT . 0 . Q5 ) . AND . ( PHD . LE . 0 . E3 ) )GO TO 401 

IFC (PHD . GT . 0 . C3 J .AND. ( PMD , LC . 0 . 30 ) ?GO TO 403 
IF < t THD . OT , 0 . 30 >. AND . t THD . LE . 1 . 00 > )GO TO 403 
IF( (PHD , GT . 1 . 00? .AND. (FHD.LC. 1 .30) 1GO TO 407 
IP C t PHD . GT . 1 . 30 > . AND . CFHD . LE . 3 . 00 ) >GO TO 408 
IPtPHD.OT, 3. 00)00 TO 411 
401 CDO* 0 .63+0 . C3*trHD~0 .05) 

403 CDO* 0 . 6B+0 . 1 C* C PHD- 0 . E5) 

403 CDO* 0 . T1 +0 . 03* t PHD"* 0 . 5 ) 

40T CDO*0 .TE3+0 . 01 ♦CPHD-l . 0 ) 

408 CDO*D .73+0 . 0D53*CPHD-1 .5) 

411 CDO* O'. T38 

FK1*£. 0/3.0 

FKE* ( £ . 0*9 .81 >**0 .5 

U* 1.5 

FU*1 . 0/U 

HDMAX*CQBP1AXXC CD0*PK1 *FKC > 1**PU 
HI MAX®HDMAX-VA 
r*TlMAX-HlMAX 
111 CONTINUE 

HA*CVA**£. 0)/C£. 0*9.81 ) 

H1MIN=CQBMIN/C CD0*FK1 *FK£ ) )**FU 

P*Y1MAX-H1MAX 

WRITER 6, 43) 

48 FORMAT U OX, * DESIGN PARAMETER ') 

WRITE (6, 1 04 )P> HI MAX, HI MIN 

t 04 FORMAT < 1 OX , ' HEIGHT OF CREST ABOVE APPROACH CHANNEL BED P'/ 

1 1 OX, ' IN METERS'/ 

1 t OX, 'DESIGN HEAD ABOVE CREST OF WEIR H1MAXV 
1 1 OX , ' MINIMUM HEAD ABOVE CREST OF WEIR HI MIN * / 

1 10X/ P« ' FI 0 .4/1 OX, 'HI MAX* 'FI 0 .4/1 OX, ' HI MIN* ' F 1 0 . 4/ ) 

I F ( IANS .EG . 1 ) GO TO 331 
GO TO 33E 

331 WRIT£(E0,38)QinAX,QlMIN,B, Y1MAX 

WR I TE ( £ 0 i 48 ) 



u u 


UR I TE i £0 , f 04 )F,H!MAX,HtMIN 

UR I TE C 2 0 , * ) ' OUTPUT FOLLOWS- — 

332 WRITEC6,*)' OUTPUT FOLLOWS # 

X = 0 . 00 

DX = 0 . 1 

502 XH*»X/HDP1 AX 

YD«HDMAX* 0 . 5* t XH ) * * 1 .25 
WRITEC6, 566 ) X, TD 
IFC IANS . EG. 1 )GO TO 333 
GO TO 334 

333 WFIT2<20, 566) X, TD 

566 FORMAT C 1 OX, ' X« ' , F 1 0 . 5 , 2X , ' YD= ' , FI 0 . 5 ) 

334 X*X+DX 
AS=r 

XFCTD.GT.A5) GO TO 501 
GO TO 302 
301 X=. 00 

DX=0 . 1 

504 XH=X/HDMAX 

TU=HDMAX* ( 0 . 724* < XH+C . 2T ) ** 1 .25- 
1 . 432* C XH + . 2T ) **0 . 623+ .126) 

WKIT2C6, 3 6 T ) X, TU 
IF( IANS . EQ . 1 >GO TO 335 
GO TO 336 

333 WRIT2C20, 36T) X, TU 

56T FORMAT C 1 OX , # X= ' , F 1 0 . 5 , 2X , ' TU= ' , FI 0 . 5) 

336 X=X+DX 

25=0 .270*HDMAX 
IFC TU . GT . 22 )GO TO 503 
GO TO 504 
503 RETURN 

END 


ROUND NOSE HORIZONTAL BROAD CRESTED WEIR 
SUBROUTINE RNBW(IANS) 

WRITEC6, 1 00) 

100 FORMAT ( t OX, 'WELCOME TO THE DESIGN OF '/ 

1 1 OX., ' ROUND NOSE HORIZONTAL BROAD CRESTED WEIR'/) 

WRITE* 6,32) 

32 FORMAT C 1 OX , * ENTER Qt MAX, Q1 MIN, B, FL, Y1 MAX ' / ) 
READ<6,*)QtMAX,Q1MIN,B,FL,YlMAX 

PAUSE 

UR I TEC 6, 33) 

33 FORMAT* 1 OX, 'ENTER TRIAL VALUE OF Ht MAX ' / ) 

READ ( & , * ) HI MAX 

WRITE (6 ,38)Q1MAX,GtMlN,B,FL, Yt MAX 


38 FORMAT ( t OX , 1 INPUT DATA V 

1 t OX , ' MAXIMUM DISCHARGE Q1 MAX= ' , FI 0 . 4 , 2X , ' CUBIC METERS 9 / 



I 1 V ^ f nininvil 74** * mwm wrrn* w **»«»«*. « - - ... . „ .... 

1 10X, 'BREATH OF WEIR B«',FtO 4,£X, 'METERS ' / 

t 1 OX, 'MAXIMUM DEPTH OF WATER Y1HAX*' ,FI0.4,2X, 'METERS'/) 
X-0 003 
U-1 .5 
BM"0 . 5*8 
FKt-E . 0/3.0 
FK2*(FK1 *9 . 81 )**0.5 
DO 111 I»1 , 1 0 

R-0.2*H1MAX 

CD* C 1 -£■ . 0*X*(FL-R)/BM)*{ 1 . 0-X* < FL-R )/H1 MAX ) ** l .5 
A*B* Y 1 MAX 

VH* ( 01 MAX/A )**£. 0/(2. 0*9.81 ) 

HE MX 1 *H 1 MAX+VH 

CV* ( HE MX 1/H1MAX)**U 

HI MAX* ( G1 MAX/( CD*CV*FK ? *FK£*BM ) >**0 .6666 
111 CONTINUE 

WRITE(6,48> 

48 FORMAT? 1 OX,' DESIGN PARAMETER *) 

WR I TE < 6 , 1 04)P,H1MAX 

104 FORMAT (1 OX, 'HEIGHT OF CREST ABOVE APPROACH CHANNEL BED P'/ 

1 1 OX, ' IN METERS '/ 

1 ? OX, 'DESIGN HEAD ABOVE CREST OF WEIR HI WAX ' / 

1 1 OX, 'MINIMUM HEAD ABOVE CREST OF WEIR HI MIN ' / 

1 10X,' P='F1O.4/10X, 'HI MAX* ' F 1 0 . 4 / 1 OX/ ) 

P = Y1 MAX-H1 MAX 
H 1 M I N= 0 .E 
DO Till 1=1,5 
R=0 . E*H? MIN 

CD* ( 1 -E . 0*X* (FL-R )/BM)* ( 1 . 0-X* ( PL-R ) /HI MIN ) ** 1 . 5 
Q 1 MI N® 7 .0 

H1MIN*( 0?MIN/( CD*rKl *rKE*BM) 7**0.66666 
1111 CONTINUE 

WRITE ( 6 , 35 ) H1MIN 
35 FORMAT! 10X, ' H 1 MI N= ' , Fl 0 . 4 ) 

WRITE (6, 41) 

41 FORMAT! 1 OX, 'LIMITATION OF ROUND NOSE BROAD CRESTED WEIR*/ 

1 1 OX, 'HI MAX? 0 . 06 ' /I OX, 'HI MAX>0 . 05*FL'/ 

1 1 OX , ' P > 0 . 1 5 METER'/lOX, ' H 1 MAX/P< 3 . 0 ' / 

1 1 OX, ' 0 . 05< ( H1MAX/FL) >0 . 5'/l OX, 'B?0 .3'/ 

1 10X, ' B? HI MAX ' / 1 OX , 'B>0.2*FL'/J 

I F ( IANS.EQ. 1 ) GO TO 339 
GO TO 340 

339 WRITEtEO, 38)01 MAX ,Q1MIN,B,FL,Y1MAX 
WRITEC20, 4S) 

WR I TE ( 6 , ? 04 ) P , H 1 MAX 
WRITE (20,35) HI MIN 
WRITE (6,4? ) 

340 RETURN 
END 

C 

C ROMJIN MOVABLE WEIR 

SUBROUTINE RMU(IANS) 

WRITE (6,100) 

100 FORMAT ( 1 OX , ' WELCOME TO THE DESIGN OF'/ 

1 1 0 X , ' ROM J I N MOVABLE WEIR '/) 

PAUSE 

WRITE(6, 1032) 

1 032 FORMAT? 1 OX, ' Q1 MAX , Q 1 MI N , Y1 MAX, FL , B ' / ) * 



n n 


READ(5, *>Q1MAX,CMMIN, Y1MAX,FL,B 
PAUSE 

WP I T E C 6 , 36 ) O I WAX , Q J « I N , B , T IP AX 
38 FORMAT C 1 OX , ' IN PUT DATA '/ 

1 1 OX , ' MAXIMUM DISCHARGE Q 1 MAX* ' ,F10.4,EX, ' CUBIC METERS'/ 

1 1 OX, 'MINIMUM DISCHARGE Q 1 M I N= ' , P 1 0 . 4 , EX , 'CUBIC METERS'/ 

1 10X, 'BREATH OF WEIR B* ' , FI D . 4, EX, 'METERS'/ 

J 1 OX, 'MAXIMUM DEPTH OP WATER Y 1 MAX* ' , P 1 0 . 4 , EX , ' METERS ' / ) 

WRITEt 6, 1 036 ) 

1036 FORMAT ( I OX , ' ENTER TRIAL VALUE OF H 1 MAX ' ) 

READ C5,*)H1MAX 
DO 111 1=1,10 

CD= 1 . 0 
A=B* Y 1 MAX 

VH=t Q1 MAX/A )**£.0/C£.0*?.Bl ) 

HEMX1 =H1 MAX+VH 
0=1.5 

CV= CHEMX1 /HI MAX ) **U 

FK 1 =E . 0/3 . 0 

FK£= 1 PK1 *9 . SI ) **0 . 5 

HI flAX x t CM MAX / C CD*CV*PK1 »PKE*BM J J **0 . 6666 
HtniN=tail1XN/(CD»CV*rK1*FKE*Bn) )**0.6666 

111 'continue 

WR I TE( 6 , 4B ) 

WRITE(E0,4S> 

A8 FORMAT ( 1 OX , ’ DESIGN PARAMETER *) 

WRITE* 6, 1 04 >P, HI MAX, HI MIN 

10A FORMAT* 1 OX, 'HEIGHT OF CREST ABOVE APPROACH CHANNEL BED P'/ 

1 1 OX, ' IN METERS'/ 

I 1 OX, 'MAXIMUM HEAD ABOVE CREST OF WEIR H1MIN IN METERS'/ 

1 10X, 'MINIMUM HEAD ABOVE CREST OF WEIR HI MIN IN METERS'/ 

1 10X,' F= 'FT 0 . 4/1 OX, 'HI MAX*'F1 0 .4/1 OX, 'H1MIN*'F1 0 .4/) 

342 RETURN 

END 

*%**************** ********************************************4 

TRIANGULAR BROAD CRESTED WEIR 
SUBROUTINE TBW(IANS) 

WRITE (6,100) 

100 FORMATC 1 OX, 'WELCOME TO THE DESIGN OF'/ 

1 1 OX, 'TRIANGULAR BROAD CRESTED WEIR'/) 

PAUSE 

WR I TE (6,11 3£ ) 

113£ FORMATUOX, ' Q1 MAX, Q1 MIN, Y1 MAX, B, ANG, FL'/) 

READ ( 6 , * ) Q 1 MAX , Q 1 M I N , Y 1 MAX , B , ANG , FL 
WR I TE (6,1136) 


FORMAT* 1 OX, 'ENTER TRIAL VALUE OF HI MAX ' ) 

READ ( 6 , * )H1 MAX 
U = £ . 5 
P=1 . 0 
FU= 1 . 0/U 
A=B* Y 1 MAX 


t 136 



UU ! I I 1*1/5 

AL*H 1 MAX /FL 

IT ( (AL C7 . 0 . 06 ) . AND . ( AL . LE . 0 . 1 > ) GO TO 410 
IFUAL CT . 0 . 1 0 ) . AND . ( AL . LE . 0 . 1 5 ) > CO TO 412 
ire ( AL GT. 0.15). AND. (AL.LE.O.£))GO TO 414 
IF ! ( AL CT.O.SO) .AND. ( AL . LE . 0 . 3 ) ) CO TO 416 
IF t AL . GT . 0 . 30 ) CO TO 418 
410 AHL*AL- 0.06 

CD* 0 . 85+1 . 30 + AHL 
4 1 E AHL*AL-0 . 1 0 

CD* 0 .902 + 0 . 64+AHL 
414 AHL*AL-0 . 15 

CD*0 . 93E+0 . 32+AHL 
416 AHL-AL-0.2 

CD=0 .948+0 . 15+AHL 
418 AHL*AL-0.3 

CD=0 . 963+0 . 1 *AHL 

HB=0 .5*BM 

HEMX1=H1MAX+VH 

WR I TE !*,*)' HEMX1 ' , HEMX1 

CV- ( HEMX1 /HI MAX ) * *U 

FK 1 * 0 . 800 

FK£ = ( 0 . 4 + 9 . 81 )**0 . 5 
THETA=(3. 14*ANG)/180. 0 
FK3=TAN! THETA ) 

HI MAX* < 01 MAX/ ( CD*CV*FK1 +FKE+FK3 ) )**FU 
HI MIN* ( C31 MINX ( CD+OV+TK1 +FKE+FK3 ) ) **FU 
34 FORMAT! 1 OX , 'H1MAX* ' , FI 0 . 3) 

111 CONTINUE 

F=T1MAX-H1MAX 
Wft I TE t 6 , 48 ) 

WRITE CEO, 48) 

48 FORMAT! 1 OX,* DESIGN PARAMETER ’) 

WRITE* 6, 1 04 >P, HI MAX, HI MIN 

104 FORMAT! 1 OX, 'HEIGHT OF CREST ABOVE APPROACH CHANNEL BED F'/ 

1 1 OX, 'IN METERS'/ 

1 IOX, 'DESIGN HEAD ABOVE CREST OF WEIR H I MAX ' / 

1 IOX, 'MINIMUM HEAD ABOVE CREST OF WEIR HI MIN ' / 

1 IOX,' P='F10 .4/10X, 'H1MAX*'F1C.4/10X, ' H 1 N I N= ' F 1 0 . 4/ ) 

IF! IANS . EQ . 1 )GO TO 343 
GO TO 344 

343 WKITEt CO, 343)01 MAX, Q1MIN,Y1MAX,B,ANG 

345 FORMAT! IOX,' INFUT DATA '/IOX, 

1 IOX, 'MAXIMUM DISCHARGE Q1 MAX* ' , F 1 0 . 4 , EX , ' CUB I C METERS'/ 

1 IOX, 'MINIMUM DISCHARGE Q I MI N= ' , F 1 0 . 4 , 2X , ' CUBI C METERS'/ 

1 IOX, 'MAXIMUM DEFTH OF WATER Y 1 MAX* ' , F 1 0 . 4 , EX, ' METERS ' / 

1 IOX, 'BREATH OF WEIR B= ', F 1 0 . 4 , £X , ' METERS '/ 1 OX ,' ANGLE* ' , F5 .E/3 
URITE(20, 48) 

WRITE! 6 , 1 04 )F , HI MAX , HI MIN 


344 


RETURN 
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C BROAD CRESTED RECTANGULAR PROFILE WEIR 

SUBROUTINE BCRCIAN5) 

WRITE (6, 100) 

100 FORMAT (I OX , ' WELCOME TO THE DESIGN OF BROAD CRESTED'/ 

1 1 OX, 'RECTANGULAR PROFILE WEIR'/) 

PAUSE 

WRITEC6, J01 J 

101 FORMATM OX, 'ENTER Q 1 MAX , Q 1 M I N , B , Y 1 MAX , FL ' , / ) 

READ ( 5 , * ) Q 1 MAX , Q I MI N , B , T 1 MAX , FL 
WRITEC6,380)Q1MAX,Q1MIN,B, Y1MAX 

380 FORMAT t l OX , IN PUT DATA •/ 

1 t OX, 'MAXIMUM DISCHARGE Q 1 MAX- ' , F 1 0 . 4 , EX , ' CUB I C METERS'/ 
1 10X, 'MINIMUM DISCHARGE Q 1 M I N = ' , F 1 0 . A , EX , ' CUB I C METERS'/ 

1 10X, 'BREATH OF WEIR B= ' , F 1 0 . 4 , EX , ' METERS ’ / 

I 10X, 'MAXIMUM DEPTH OF WATER T 1 MAX= ' , F 1 0 . 4 , EX , ' METERS ' / ) 
WR I TE ( 6 , 1 02) 

10E FORMAT* 1 OX, 'ENTER TRIAL VALUE OF H 1 MAX ' / ) 

READ* 5, * )H1 MAX 
U= 1 .5 

DO 111 1=1,10 

P=Y1 MAX-H1 MAX 
AFL*H1 MAX/FL 
APL=H1MAX/CH1MAX+P ) 

IFCAFL. LE.O. 35 ) GO TO El 

IFC CAPL.GT.0.35) .AND. (APL. LE . 0 4 0 ) ) GO TO EE 
IFC C AFL . GT . 0 . 40 ) . AND . C AFL . LE . 0 .50 ) )GO TO E4 
IFC C AFL .GT . 0 . 50 ) . AND . C AFL . LE . 0 . 55 ) )GO TO E6 
IFCAFL. GT.O. 55 ) GO TO ET 

51 irt AFL . LE . 0 . 4 )GO TO EB 
IFt AF-L . GT . 0 . 4 )GO TO 30 

SB F=1 . 0 

GO TO 50 

10 F=1 . O+C 0 . E4/1 . 1 )*C AFL-0 . 4) 

GO TO 50 

!£ IFC AFL . LE . 0 . 35 ) GO TO 3E 

IFCAFL.GT.0.35)GO TO 34 

52 F= 1 .01 
GO TO 50 

14 F = 1 . 01 + . 1 6* t AFL-0 . 35 ) 

GO TO 50 

!4 IFCAFL .LE. 0 ,35)G0 TO 36 

IFCAFL. GT.O. 35 ) GO TO 38 
;6 F = 1 . 02 

GO TO 50 

e F=1 . 02+0 . 18*CAFL-0 . 35) 

GO TO 50 

6 IFCAFL. LE.0.35)GO TO*40 

IFCAFL. GT.0.35)GO TO 42 


10 F = 1 . 04 

GO TO 50 

2 F=1 . 04+0 . 2* ( AFL—0 . 35 ) 

GO TO 50 

7 IFCAFL. LE.0.35)GO TO 44 



I3C 


46 

50 


1 1 1 


48 
1 06 
1 07 


1 04 


1 05 


346 


347 


CO TO 50 

P * 1 . 06+0 .E*!AFL-0.35> ' 


CO TO 50 
CD* 0 . 848*F 
A*B*(H1MAX+P> 

VH* (Q1 MAX/A )**E , 0 /Ce. 0*9.81 ) 

HEMX1 *H 1 MAX+VH 
CV*CHEMX1 /HI MAX ) **U 
FK1*E. 0/3.0 
FKE*CFK1*9 .81 )**0.5 
FU* 1 . 0/U 

HI MAX* ( Q1 MAX/ (CD*CV*FK1*FK£*B) ) **FU 
CONTINUE 

H1MIN=CQ1MIN/(0. 848*FK1 *FKE*B) ) **U 

HLMAX*H1 MAX/FL 

WR I TE ( * , * ) ' HLMAX* ' , HLMAX 


1 ) 


IFCHLMAX.GE.0.33)G0 TO 106 
WR I TE (6,48) 

FORMAT! 1 OX,' DESIGN PARAMETER_ . 

UR I TE ( 6 , 1 07) , , 

FORMAT! 1 0X, 'THE AIR POCKET SHOUD BE FULLY AIRATED /) 

P=Y1MAX-H1MAX 

WRITE (6,1 04 >P, HI MAX, HI MIN 

WRITECEO, 104)P,H1MAX,H1MIN _ a>i11ii ' D[ - n 0 , , 

FORMAT! 1 OX, 'HEIGHT OF CREST ABOVE APPROACH CHANNEL BED P / 

ASOVE CRE5T OF U EIR H 1 MAX I" "<£«*•' 

TOX, 'MINIMUM HEAD ABOVE CREST OF WEIR H1MIN IN ME 
10X, ' P*' FI 0.4/1 OX, 'HI MAX* 'FI 0.4/1 OX, ' HI M I N* FI 0 . 4/ ) 

WRITE! 6 , 1 05) _ , 

FORMAT! 1 OX, 'LIMIT OF APPLICATION OF BROAD CRESTED 

1 OX RECTANGULAR PROFILE WEIR'/ 

1 OX, 'HI MI N> 0 . 06 '/I OX, * H 1 MAX/ C H 1 MAX+P ) <0 • * ' 

10X, ' P > 0 . 1 5 ' / 1 OX , 'HI MAX/FL < 1.5'/10X,'BM>0.3 /) 

IF! IANS .EQ . 1 )GO TO 346 
GO TO 347 

WRITE CEO, 38 0 ) Q 1 MAX , Q 1 M IN , B , Y1MAX 
WRITECEO, 48) 

WRITECEO, 1 07) 

WRITE (SO, 1 04 )P, HI MAX, HI MIN 

WRITECEO, 1 05) 

RETURN 

END 
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EXAMPLE 1 

Q1MAX-MAXIMUM DISCHARGE IN CUBIC METERS/SEC 

G1 MIN=MINIMUM DISCHARGE IN CUBIC METERS/SEC 

B-WIDTH Or WEIR IN METERS 

HI MAX*MAXIMUM HEAD IN METERS 

rL=LENTGH ALONG TLOW DIRECTION 

YD=DOWN STREAM T COORDINATE PROFILE 

OF WES~SF ILLWAY METERS 

YU-DOWN STREAM Y. COORDINATE PROFILE 
OF WES-SPILLWAY IN METERS 

INPUT DATA WES STANDARD WEIR 

MAXIMUM DISCHARGE G1MAX* 60.0000 CUBIC METERS 

MINIMUM DISCHARGE CM M IN* 6.0000 CUBIC METERS 

BREATH OP WEIR B= 15.0000 METERS 
MAXIMUM DEPTH OF WATER Y 1 MAX* 8.0000 METERS 

DESIGN PARAMETER 

HEIGHT OF CREST ABOVE APPROACH CHANNEL BED P 
DESIGN HEAD ABOVE CREST OF WEIR HI MAX 
MINIMUM HEAD ABOVE CREST OF WEIR HI MIN 
P= 7. 0005 METERS 

H1MAX* .5551 METERS 

HI MIN* .3230 METERS 


WEIR PROFILE 


x= 

. 00000 

YD* 

. 00000 

x= 

. 1 0000 

YD* ’ 

. 00501 

x* 

. £0000 

YD* 

.01605 

x= 

. 30000 

YD= 

. 038E1 

X- 

. 40000 

YD= 

. 06506 

X- 

. 50000 

YD= 

. 09831 

x= 

.60000 

YD* 

. 13775 

x= 

.70000 

YD= 

. 183E1 

x= 

.80000 

YD* 

.£3455 

x* 

. 90000 

YD = 

.£9165 

x= 

1 . 00000 

YD= 

. 3544£ 

x= 

1 . 1 0000 

YD= 

. 4££76 

x= 

1 . £0000 

YD= 

.49659 

x* 

1 .30000 

YD* 

.57585 

x= 

r. 40000 

YD= 

. 66047 

x= 

1 .50000 

YD= 

.75038 

x= 

1 . 60000 

YD= 

.84554 

x= 

1 .70000 

YD= 

. 94590 

x= 

1 .80000 

YD* 

1 . 05140 

x= 

1 . 90000 

YD* 

1 . 1 6£0 1 

x= 

£.00000 

YD* 

1 . £7767 

X- 

£ . 1 0000 

YD* 

1 .39836 

x= 

£.£0000 

YD* 

1 . 5E404 

x= 

£.30000 

YD* 

1 . 65467 

x= 

£.40000 

YD* 

1 . 79 0£ 1 

x* 

£.50000 

YD* 

1 .93065 

x= 

£.60000 

YD* 

£.07594 
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X* 

e. 7oaoo 

YD* 

£ . ££606 

X* 

E.aoooo 

YD* 

£.38098 

x* 

£ . 90000 

YD* 

£.54068 

x* 

3 . 00000 

YD* 

£.70513 

x* 

3. 1 0000 

YD* 

£ .87431 

X* 

3. £0000 

YD* 

3 . 04819 

X * 

3 . 30000 

YD* 

3 .££675 

X* 

3 .40000 

YD* 

3.40997 

x* 

3.50000 

YD* 

3 . 5978E 

x* 

3 .60000 

YD* 

3.79030 

X * 

3 .70000 

YD* 

3 . 98738 

X = 

3 .80000 

YD* 

4. 18903 

X* 

3 .90000 

YD* 

4.395E5 

X = 

4 . 00000 

YD* 

4 . 60601 

X = 

4 . 1 0000 

YD = 

4.8E130 

X 3 

4 . £0000 

YD* 

5.04110 

X * 

4 . 30000 

YD* 

5 .£6539 

X * 

4 .40000 

YD* 

5 . 49416 

X * 

4 .50000 

YD* 

5.7E739 

X* 

4 . 60000 

YD* 

5 . 9 650 T 

x = 

4 . TOOOO 

YD* 

5 . COT1 9 

x = 

4 .80000 

YD* 

6 . 453T£ 

x* 

4 . 90000 

YD* 

8 . T 048 6 

x* 

5.00000 

YD* 

8 . 95999 

x = 

5.10000 

YD* 

T .£1969 

x* 

. 00000 

YU* 

- . 0 01 £7 

x* 

. 1 0000 

YU* 

.00506 

X * 

- . £0000 

YU* 

. 0£34 1 

x* 

- . 30000 

YU* 

, 05£85 

x* 

- . 40000 

YU* 

. 09£7 1 

x= 

- . 50000 

YU* 

. 14£5£ 

x* 

- . 60000 

TU* 

.£0158 

X s3 

- . TOOOO 

YU* 

.£7045 

x* 

- . 50000 

YU* 

. 34507 

x* 

- .90000 

YU* 

. 43443 

EXAMPLE 

£ 




cnnAX=tiAxmufi discharge in cubic pieters/sec 
G nniN*=MiNiiiurn discharge in cubic pieters/sec 

B-W IDTH OF WEIR IN PIETERS 
HIPIAXsflAXIMUn HEAD IN PIETERS 
^PROPORTIONALITY PARAMETER 

X = X COORDINATE OP PROPORTIONAL WEIR IN PIETERS 
T = Y COORDINATE OP PROPORTIONAL WEIR IN PJETERS 

INPUT DATA LINEAR PROPORTIONAL WEIR 

Q 1 MAX= . A 0 0 0 Q1MIN= .1500 B= .6000 

T1P1AX* 1.0000 CD= .6300 G= 5.B100 

DESIGN PARAMETER 

,£6£A K= .85761 H1MAX= 


A- 


.55387 P 


.4461 



Y* 

. 01 000 

ym 

. 0 £ 0 0 0 

r * 

. 03000 

T* 

. 04000 

Y= 

. 05000 

Y* 

. 06000 

Y = 

. 07000 

T r 

. 05000 

T = 

. 05000 

T= 

. 1 0000 

Y = 

. 1 1 000 

Y* 

. 1 £ 0 0 0 

T = 

. 13000 

T* 

, 14000 

T = 

. 15000 

T- 

. 16000 

T = 

. 17000 

Y = 

. 15000 

Y = 

. 1 5000 

Y* 

. £0000 

Y*= 

. £1 000 

Y = 

. ££000 

T* 

. £3000 

Y E 

. £4000 

7 s 

. £5000 

Y c 

. £6000 

Y«= 

. £7000 

Y«= 

. £5000 

Y®= 

. £5000 

Y = 

. 30000 

Y s 

.31000 

Y*= 

. 3 £ 0 0 0 

Y s 

. 33000 

Y = 

. 34000 

Y*= 

. 35000 

Y* 

. 36000 

Y*= 

.37000 

Y = 

. 35000 

Y = 

. 35000 

Y« 

.40000 

Y = 

.41 000 

Y= 

. 4£ 0 0 0 

Y* 

. 43000 

Y = 

. 44000 

Y= 

. 45000 

Y= 

. 46000 

Y* 

.47000 

Y = 

.48000 

Y== 

.49000 

Y~ 

.50000 

Y = 

.51000 

Y = 

. 5£00Q 

Y= 

.53000 

Y= 

. 54000 

Y« 

.55000 


. £6 3 1 6 
. £4555 
. £ 377 £ 
.££550 
.££135 
. £1451 
. £0554 
. £0364 
. 1 5551 
. 1 5437 
. 1 50£5 
. 1 5643 
. 1 5£36 
. 17551 
. 1 7635 
. 17335 
. 17056 
. 1 6755 
. 1 6534 
. 1 6 £ 5 £ 

. 16061 
. 15540 
. 1 56£5 
. 154£5 
. 1 5£3 0 
. 15043 
. 14563 
. 14655 
. 1 45 ££ 

. T 436 0 
. 1 4£04 
. 14053 
. 13507 
• 13766 
. 1 36£5 
. 13456 
. 13367 
. 1 3 £ 4 1 
.13115 
. 13001 
. 1 £886 
. 1 £77 4 
. 1 £665 
. 1 £558 
. 1 £455 
. 1 £354 
. 1 ££55 
. 1 £ 1 55 
. 1 £065 
. 11573 
. 11883 
. 11755 
.11710 
. 1 1 6£6 
.11544 


WEIR PROFILE 
X * 

X* 

X = 

x*= 
x~ 
x= 
x= 

X s 

x= 
x= 
x= 
x= 
x= 
x= 
x= 

X- 

x= 
x= 
x= 
x= 
x= 
x= 
x= 

X s 

x= 
x= 
x= 

X- 

x= 
x= 
x- 
x* 
x= 
x= 
x= 
x= 
x= 
x~ 

X- 

x= 
x= 
x= 
x= 

X- 

x= 
x= 
x= 
x= 
x~ 
x= 
x= 
x= 
x= 
x= 


EXAMPLE 3 

Q1HAX*MAXIMUM DISCHARGE IN CUBIC METERS/SEC 
Q1MIN*MINIMUM DISCHARGE IN CUBIC METERS/SEC 
B*UIDTH OF WEIR IN METERS 
H1MAX-MAXIMUM HEAD IN METERS 
K-PROPORTIONALITY PARAMETER 

X = X COORDINATE OF PROPORTIONAL WEIR IN METERS 
Y*Y COORDINATE OF PROPORTIONAL WEIR IN METERS 

INPUT DATA LOGRITHMIC WEIR 

Q 1 MAX* .1100 Q1MIN* .0500 B* .5000 

Y 1 MAX* .4000 CD* .6300 G* * 9.6100 

DESIGN PARAMETER 

A* . £395 K= .08250 HI MAX® .52117 P- -.1212 


WEIR PROFILE 


Y= 

. 05000 

X = 

. 09823 

Y = 

. 1 0000 

x= 

. 06003 

Y = 

. 15000 

x= 

. 04002 

Y* 

. £0000 

x= 

. 02798 

Y*= 

. £5000 

x*= 

. 0201 6 

T* 

. 30000 

x~ 

. 01 481 

Y = 

. 33000 

x= 

. 01 1 01 

T® 

. <90000 

x~ 

. 00823 

T* 

. *33 0 0 0 

x« 

. 006 1 4 

T® 

.30000 

X s 

. 00436 

Y = 

. 33000 

x= 

. 00333 

Y« 

.50000 

X*® 

. 00236 

T s 

. 53000 

x= 

.00139 

Y* 

. TOOOO 

X s ® 

. 00098 

Y = 

. T5000 

X s ® 

. 00049 

Y* 

. 60000 

x- 

. 00008 



